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ABSTRACT 



The objective of this study was to describe the mechanics of wind wave 
generation and spectral development. Intermittency , high frequency micro- 
structure in wind and wave fields, and strong nonlinear coupling involving 
a wide range of scales are shown to be crucial elements in the transfer of 
momentum to, from, and within the wave field. None of these elements are 
included in available theories. 

Measurements of wave height and of the turbulent atmospheric and sub- 
surface boundary layers were made from a small surface following platform 
and from a stable 38.5m spar buoy. The structure of moving gust patterns 
(cat’s paws) is described and related to the generation of surface waves. 
Results from this and other background studies are then applied to a dis- 
cussion of spectral growth during a two day period of active wave 
generation. 

Cat's paws contain 'bursts' of intense turbulent stress and buoyancy 
fluctuations separated by quiescent 'intervals'. There is a difference of 
over three orders of magnitude in fluctuation strength between these 
features. Rapid growth rate generation of high frequency surface waves 
and atmospheric turbulence occurs during the bursts. The resultant micro- 
scale components aid the growth of lower frequency instabilities by strong 
nonlinear coupling between scales of motion and by acting as drag or 
roughness elements. Evidence of strong coupling between frequency bands 
and of weakly resonant capillary-gravity wave interactions is presented. 
Thermal stratification has a strong influence on fluctuation magnitude 
and can delay the onset of surface wave generation. 

Major spectral growth is highly unsteady. Much of the momentum flux 
from air to sea occurs during intermittent events that are similar in 
nature to cat's paws, and goes directly into high frequency waves. The 
bursts occur predominantly over large groups of surface waves and involve 
strong nonlinear interactions between media and frequency bands. The 
long-term equilibrium balance between wind and water is disrupted by 
variations in surface currents. There are 'critical' wind speeds charac- 
terized by anomalous relationships between parameters of the predominantly 
logarithmic velocity profile. 
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I. INTRODUCTION 



Description of the Study 

This paper describes an experimental study of the mechanics of wind- 
wave generation and the spectral development of sea surface waves. Em- 
phasis is placed on the time history of generation and development, and on 
the nature and time scales of air-sea coupling. The basic objective was 
to describe and illustrate the physics involved in the transfer of 
momentum to, from, and within the wave field. The underlying hypothesis 
is that this process is highly nonlinear and depends strongly on the 
microstructure of the wave and wind fields. The experiments were con- 
sequently focused upon investigation of intermittent interactions between 
variables and scales of motion. 

Two separate experiments are described. The first is a detailed 
investigation of the structure of moving gusts. The second is an anal- 
ysis of spectral growth during a two day period of increasing wind. 

Both separately and together their results indicate that wind-wave gen- 
eration and spectral development involve strong nonlinear interactions 
with the wind and current fields; representation of these processes in 
terms of linear or weakly nonlinear coupling is inadequate. 

Background 

The available linear and weakly nonlinear wave generation models 
fail to yield satisfactory estimates of wave growth. In addition there 
exist unexplained anomalies such as spectral overshoot and the degener- 
ative effect of surface slicks on long waves. This combination of qual- 
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itative and quantitative unpredictability indicates that the physics of 
wave generation and spectral development is still inadequately under- 
stood. Until the correct physics is described, there can be little sig- 
nificant progress in either theoretical modeling or forecasting. 

Since 1957, most wind-wave generation experiments have been designed 
around the linear Miles-Phillips models (Miles [1957], Phillips [1957, 
1966]). Early studies compared observed spectral growth with theoretical 
values; they soon showed that spectral shape and initial growth rates 
were reasonably well predicted, but that growth rates of the major por- 
tion of the spectrum were severely underestimated. Consequently, more 
recent work in the low frequency range has concentrated on the direct 
relationship between wind and wave fields in attempts to determine the 
adequacy of the theory’s representation of momentum transfer mechan- 
isms. The high frequency range (f > lhz) has been treated as a separate 
entity, primarily relevant to initial generation. 

Work based on the simplified linear model has essentially dealt 
with the general or time- averaged characteristics of wind-wave coupling 
and spectral development. It has neglected the relationship between 
portions of the wave spectrum, and in particular has deemphasized the role 
of high frequency wavelets and the coupling between waves and turbulence. 
There has thus been very little work encompassing the full range of spec- 
tral growth and relating it to interactions within and among the wind, 
wave, and surface current fields. Hasselmann * s Hamiltonian type model 
(1967) was so general (and unfamiliar to field workers) that it also 
failed to inspire experimental investigation along these lines (with 
the exception of some ongoing parts of JONSWAP) . 

A second experimental approach to the problem of spectral growth 
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has been through parameterization of wind and wave fields and the devel- 
opment of spectral transfer functions for use in forecasting. Such studies 
have described the generally logarithmic nature of the atmospheric 
boundary layer over water, specified the characteristics of wave spectra 
in some detail, indicated the significance of energy transfer between 
wave components, hinted at the importance of high frequency elements, 
and yielded partially valid prediction schemes. They provide little 
direct information about the physics of momentum transfer processes, but 
rely upon theory and experience for guidance in the choice of empirical- 
ly applied parameters. 

It is highly unlikely that either forseeable modifications to the 
available models or additional parametric type studies will correct the 
known theoretical deficiences. Further advances will require a more 
complete understanding of the physics of generation. Important elements 
in a reevaluation of the wave generation process must include the detailed 
structure of the turbulent boundary layers and nonlinear couplings, both 
weak and strong, within and among the three fields of motion. Theoreti- 
cal and experimental studies have shown that such couplings do exist. 
Turbulence is known to be intermittent and strongly dependent on boun- 
dary conditions; a large number of scales are involved in the dynamics of 
flow at typical geophysical Reynolds numbers; and the generation of waves 
is inseparable from the generation of turbulence. These considerations 
imply that wave and turbulence generation may be a strongly nonlinear pro- 
cess involving a large variety of forces and scales, and that the evalu- 
ation of experimental results in terms of mean values and long averaging 
times can easily obscure the physics. Viewed in these terms the time his- 
tory of events during wave generation becomes an important factor; the 
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microstructure of the air and water boundary layers and the wave field 
reassumes great significance; and the ability to distinguish the differ- 
ences between generation, transfer, and detritus becomes a major experi- 
mental goal. 

The pair of experiments described herein was designed to investigate 
the waves and near-naviface boundary layers and to determine the nature 
and time scales of the important processes in wave generation and spectral 
development. Events ranging in duration from less than one second to 
two days were examined to obtain answers to the following questions: 

1. What are the important time scales in the transfer of 
momentum to, from and among surface waves? 

2. What are the nature and magnitude of the fluctuations 
involved in these transfers? 

3. What is the nature of the coupling between media, and 
between scales of motion both within and between media? 

4. What is the ’imprint 1 of wave generation in the air and 
in the water; and can this imprint be used to distinguish 
events of importance in the generation process (i.e., both 
to predict and to permit conditional sampling)? 

5. Is there a similarity between scales, and if so, for what 
variables and over what range of scales does the similarity 
exist? 

The first experiment was designed to isolate and intensively analyze 
an event of known significance - the ’cat’s paw*. The cat’s paw is an 
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isolated parcel of atmospheric turbulence. It contains anomalously large 
velocity and temperature fluctuations, and can be visually delineated by 
the 'dark 1 , moving patch of high frequency waves which it creates. The 
second experiment yielded data from a larger system. It was designed to 
extend the time scales and number of forces involved to values typical 
of open ocean spectral development. The experiments are first individually 
discussed in separate papers; then the above questions are answered and 
the development of a wave field is described in terms of the results. An 
outline of experimental details is contained in appendices; complete des- 
criptions and all significant data are available in additional reports of 
the M.I.T. Buoy Project. 
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II. THE STRUCTURE OF MOVING GUST PATTERNS (CAT T S PAWS) 
AND THEIR ROLE IN WAVE GENERATION 



Abstract 



Measurements of wave height, turbulent velocity, and temperature 
fluctuations within 1 meter of the sea surface were taken from a small 
surface following platform. The resultant data show that moving gusts 
are highly structured, consisting of bursts of turbulent stress separ- 
ated by relatively quiescent intervals. Magnitudes of stress and buoy- 
ancy fluctuations vary between bursts and intervals by over 3 orders of 
magnitude. The bursts are associated with the generation of very high 
frequency wavelets which in turn serve to couple the wave field .to the 
mean air flow. A large range of scales of motion interact, leading to 
a lognormal distribution of turbulent velocity fluctuations and to ener- 
gy transfer among wave components. 



Introduction 

Gust patterns moving over the sea surface (cat’s paws) are well- 
known to sailors, both those seeking a strong wind and those trying not 
to capsize. Cat’s paws can be located by their imprint on the water 
surface, namely the change in surf ace reflectivity caused by the high fre- 
quency wavelets generated. Associated with these gusts are large magni- 
tude, high frequency, turbulent stress and buoyancy fluctuations. Cat’s 
paws are patches of intense air-sea Interaction with correspondingly large 
transfers of momentum, heat, and moisture. This study reports an attempt 
to measure, describe, and find relationships among the fluctuations with- 
in a cat’s paw, with particular reference to the mechanisms involved in 
momentum transfer and wave growth. 

Incorporated in a cat’s paw are three elements of wave generation 
and spectral development that are neglected in the available modern models 
(Jeffreys [1925, 26], Phillips [1957, 1966], Miles [.1957], Hasselmann 
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[1967]). These are intermittency of generation, the role of wavelets in 
spectral development, and direct, strong nonlinear coupling between scales 
of motion and between the sea surface and atmospheric turbulence. The cat’s 
paw is itself an intermittent or burst type phenomenon; the activity with- 
in it is qualitatively and quantitatively different than that observed 
on the average. As an experimental objective, a gust serves as a realiza- 
tion one type of wave generation event, with constraints -determined by the 
environment, at Reynolds numbers unattainable in the laboratory. 

Although the above three elements have been largely neglected, sub- 
stantial eveidence of their significance is available in the literature. 

Munk noted in 1955 that while high frequency waves account for an insig- 
nificant portion of total wave energy, they play a disporportionate role 
in both wind stress and the development of the spectrum. This had been 
exemplified by Van Dorn’s (1953) experiment, where the presence of a sur- 
face slick significantly reduced setup in an artificial pond. This 
effect of slicks was recently reemphasized (Ruggles [1969], Barger, et al, 
[1970]) by a demonstration that removal of wavelets resulted in a re- 
duction of wave energy accompanied by an increase in mean wind, decrease 
of roughness length, and significant reduction of coherence between 
wind and waves. 

Stewart (1961, 1967) has long maintained that much if not most of 
the momentum transfer from air to water is effected by wave drag. Van 
Dorn (1953) and Wu (1969 a, 1969 b, 1970) attribute this wave or form 
drag to the high frequency wavelets superimposed on larger waves. Wu 
further relates the wavelets to roughness height and states that the mo- 
mentum flux to wavelets accounts for a substantial fraction of the wind 
stress; the stress appears as a normal component to the wavelets, but 
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feels like a tangential stress to the underlying waves. Longuet-Higgins 
(1969a) has shown that such a fluctuating tangential stress, occurring 
selectively near the crests and leading faces of longer waves where the 
shorter ones are ’trapped’ (Munk [1955], Cox [1958]), is dynamically equi- 
valent to a normal pressure fluctuation in quadrature with them. 

The role of turbulence and turbulent Reynolds stress in wave gen- 
eration is commonly viewed as minor or negligible. Miles-Phillips theory 
holds that the waves do not affect turbulence; Hasselmann includes a weak 
interaction through a disturbance of mean flow by one wave component, 
interaction of the disturbance with a turbulence component, and a feed- 
back of the resulting pressure perturbation into another wave element. 
Stewart (1967) has shown that wave energy is typically comparable with 
the energy of turbulence in the entire troposphere, and consequently has 
suggested that whxle turbulence can serve as a catalyst for transfer of 
momentum from the mean flow, it can by itself add very little to the wave 
energy. 

Two of the more important properties of atmospheric turbulence with 
respect to wave generation are intermittency and the associated large 
growth rates of localized high frequency instabilities (Greenspan and 
Benney [1963]). Laboratory experiments (e.g., Klebanoff, et al; [1962], 
Kline, et al, [1967]) have shown that the progression of events leading 
to instability and transition of laminar boundary layers involves suc- 
cessively the formation of travelling Tollmien-Schlichting waves , develop- 
ment of spanwise disturbances leading to formation of local shear layers, 
and explosive development of the flow into short turbulent length scales 
in bursts of high frequency instability fluctuations. The turbulent 
boundary layer is itself intermittent, with periods of intense activity 
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separated by relatively quiescent intervals (e.g., Frenkiel & Klebanoff 
[1971], Rao and Narayanan [1971]). Reverse transition is similarly 
abrupt (Steiner [1971]) and is due to the restoring action of buoyancy 
forces. Furthermore, Klebanoff, et al, (1962) demonstrated that the 
spatial structure of turbulent breakdown is strorgxy dependent on boundary 
irregularities, and Gupta, et al, (1968) showed that unstable oscilla- 
tions iidiced by the simultaneous introduction of Tollmien-Schlichting 
waves and free surface waves resulted in a large, unpredicted amplification 
rate of small wavelengths. 

The combined effects of intermit tency , high growth rates, and boun- 
dary influence were illustrated (Mollo-Christensen [1970]) by filtering, 
squaring, and plotting selected frequency components from wave gauge and 
anemometer signals recorded during a period of wave generation. The plots 
showed that ripples and high turbulent wind energy were associated with 
groups of large amplitude, low frequency waves. This suggests that 
generation of waves and turbulence are coupled, intermittent, and asso- 
ciated with the presence of wavelets on longer wave crests. The basic 
roles of the high frequency waves are the coupling of wind and wave 
fields and the enhancement of low frequency instabilities; the wavelets 
are generated by bursts of high frequency turbulence, and growth rates 
in these local instabilities are extremely rapid and due to strong non- 
linear coupling. In turn the turbulence patterns are controlled by the 
major wave field. The combined result is a strongly coupled system with 
growth rates much larger than those predicted by linear or weak coupling 
theory. 

Implicit in this outline of spectral development is the participa- 
tion of a large range of scales, coupled by nonlinear interactions. The 
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high frequency waves are in the capillary-gravity and capillary wave 
ranges. As pointed out by Stewart and Manton (1971), the important 
time scale for generation by turbulent pressure fluctuations (Phillips* 
resonance mechanism) is that of the wave group generated, i.e., the wave 
frequency. Similarly, the large amplification rates in shear layer in- 
stability (Greenspan & Benney , [1963] ) have a typical time scale of a 
half period of the primary wave. Water surface waves of the capillary- 
gravity type have second order resonant interactions (McGoldrick [1965], 
Simmons [1969]) with effective time scales of two to three waves periods. 
Turbulent bursts or spots are similarly of short duration. It can 
therefore be anticipated that the generation of wavelets is a short- 
lived but intensive phenomenon which is capable of rapidly depleting the 
locally available energy of the mean flow. At the same time, as pointed 
out by Mollo-Chri r utensen (1971), **work done on the water to make waves must 
equal a stress times a strain, and the reactions of the stresses on the 
water act on the air to produce turbulence. *' 

The intermittency and selectivity of the location of high frequency 
waves on the underlying long ones is only partially attributable to the 
coupling with atmospheric turbulence.* Variations in surface tension near 
the crests of steep waves produce a local increase in normal stress which 
in turn generates capillary waves (Longuet-Higgins , [1963] , Crapper [1970]). 
The short waves steepen on the forward faces of longer waves due to hor- 
izontal contraction; this leads to breaking and the transfer of mass and 
momentum from the short to the longer waves (Longuet-Higgins [1969b]). 
Hasselmann (1971) has shown that the work done on the long waves by the 
interaction stresses* 1 is balanced to first order by the loss of poten- 



tial energy involved in mass transfer, and to second order results in 
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a slight attenuation of the longer waves. The basic physical idea sug- 
gested by Longuet-Higgins has not been negated by Hasselmann’s rescaling 
of the terms; there does exist a considerable energy interchange between 
scales of motion within the water, and this interchange is strongly depen- 
dent on the presence and location of high frequency waves. 

A similar interplay occurs in the atmosphere. Mollo-Christensen 
(1970, 1971) has shown that when the velocity field is described in terms 
of three components - mean flow, long wavelength and short wavelength 
fluctuations - the resultant energy equations involve interaction or pro- 
duction terms, in the form of stress times strain, among all three com- 
ponents. If the interactions occur with the proper phase relationship, 
a primary instability can be enhanced by the work done by Reynolds 1 
stresses of the secondary instability, and both instabilities can for 
a while gain energy from the mean flow. This type of interaction can 
occur similataneously on many scales, in both the wind and the wave 
fields. The essential points are that Reynolds stresses are inseparable 
from the generation of fluctuations, and that by nonlinear couplings between 
scales, the growth of short waves .can enhance the growth of longer ones. 

On the basis of these considerations from the literature, it is 
suggested that the major reason for the inability of available theories 
to satisfactorily predict spectral growth is that they fail to include 
the correct physics of the momentum transfer process. The inclusion of 
strong nonlinear couplings tremendously increases the complexity of the 
problem ( and makes it as yet mathematically intractable) , but appears 
to be necessary for a proper understanding of wave generation. 

The cat’s paw is a readily accessible phenomenon that provides- a 
realistic geophysical framework within which to initially investigate these 
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ideas. The wind tunnel is too small for studies at the Reynolds num- 
bers and hence the range of scales involved in the proposed mechanisms, 
and the dynamic range of wave height on the open ocean makes detailed 
near-naviface measurements exceedingly difficult. This experiment was 
designed to investigate the time history of cat’s paws and to determine 
the scales involved in them. The primary features sought were the rel- 
ative magnitudes of the fluctuations in air and water and an indication 
of coupling between media. Thus the variables of major concern were 
turbulent Reynolds stress and buoyancy fluctuations, and high frequency 
wave activity. Interactions between scales of motion were examined in 
terms of the basic cat’s paw structure, turbulent velocity fluctuations, 
and nonlinear wave interactions. Specifically, answers were sought to the 
following questions: 

Is there a demonstrable relationship between the generation 
of turbulence and the generation of high frequency waves? Is the 
process intermittent or is it typically near the mean or average 
flux? 



What is the evidence for coupling between scales of motion? 



On what time scales does instability-associated high growth 
rate generation occur? 



What is the imprint of a generation ’event’ in the air and 



water? 
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The Experiment 

The measurements were taken from a small (1 m long) surface following 
catamaran (cat; Figure 1). The cat was free in pitch, roll, and heave 
(design natural frequency 1.5 hz) and restrained in yaw and horizontal 
position. It could be oriented into the wind through the full range of 
acceptable fetch directions. The experiment was conducted in March 1971 
from the dock of the M.I.T. Yacht Club on the Charles River, Cambridge, 
Massachusetts. 

Instruments included 6, DISA type 55A81 hot film anemometers, a capac- 
itance wave gauge with .076 cm diameter sensor, and VECO 42A70 (air) 
and Fenwall K496 (water) thermistors balanced against precision resis- 
tors in a differential bridge circuit. Subsidiary equipment included 
an accelerometer for sensing catamaran motion and cup and propeller anem- 
ometeis for determining mean wind. Data were recorded on magnetic tape 
in analog mode, accompanied by a voice narration that was later used 
for subjective conditional digitization. 

Two of the hot films were orthogonally mounted to permit determina- 
tion of horizontal ( \J ) and vertical ( W ) velocities and hence Reynolds 
stress (-/ ? o‘ w' ). The sensors were attached to the forward vertical cat 
support to minimize tilt (<1° mean, ^3° fluctuating). Table 1 lists the 
instrumental arrangement and mean conditions for the seven 2-hour runs 
that were analyzed. Data from runs 2, 4, and 6b, which represent the 
full range of conditions encountered, were selected for presentation 
herein. 

The data were analyzed using a combination of digital and analog 
methods. Pre-analysis, demodulation, and filtering (100 hz low pass) 
were analog. The data were then digitized at 250 hz and recorded on 7 
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Figure 1 



Catamaran and Support Structure 



TABLE 1: Catamaran Runs (2 hour duration) 
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track digital magnetic tape. Final analysis was accomplished using lo- 
cally generated software on the M.I.T. IBM 360. All available channels 
were analyzed simultaneously for each run segment of interest. Exten- 
sive use was made of graphical output from a Stromberg-Carlson 4020 
plotter. 

The major contraints on the analysis were the cat motion and the 
nature of the cat’s paw itself. Cat motion could be successfully removed 
from the data, but the very fact that it did follow the low frequency 
waves prevented meaningful quantitative estimates of spectral growth. 

The cat’s paw will be shown to consist of patches of strong fluctuations 
separated by relatively quiescent intervals. Thus time averaging and 
the selection of data segments for analysis had to be very carefully 
done to avoid distortion of fluctuation magnitudes and to prevent ob- 
scuring events of importance. The time history of events proved to be 
of crucial importance. Further, the momentum transfer processes involve 
both internal and external interactions; they are highly non-s tat ionary , 
non-homogeneous , and anisotropic, so that the standard statistical theo- 
ries of turbulence are inapplicable and time averaging could not be sub- 
stituted for ensemble averaging. Therefore data from any one gust pos- 
sesses little if any statistical significance, and conclusions could be 
drawn only on the basis of the change with time in the individual prop- 
erties of a variable and in the relationship between variables, and then 
tested on the basis of similarity among a number of identically treated 
events. Each cat’s paw, and each segment within a cat’s paw, must be 
viewed as a single realization. No meaningful error bars or confidence 
limits can yet be placed on the data since a statistically insufficient 
number of events has been evaluated. The data were subjectively viewed 
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in terms of the relative changes with time, the sequence or pattern of 
events, and the relationship between variables in both time and frequency 
domains. Reported magnitudes can be considered only as order of magnitude 
estimates. This is particularly true of wave- autospectra, where filtering 
and the necessarily very short averaging times preclude even rough estimates 
of the magnitude of spectral growth and decay. 

Emphasis was placed on understanding the intermittent nature of tur- 
bulent fluctuations, their interaction with the free surface, and changes 
in the presence and structure of capillary-gravity and pure capillary 
waves. These features were then evaluated in terms of their significance 
to the generation of surface waves. 



Results and Discussion 



A. Cat’s Paw Structure in Time and Frequency . Separating velocity 
and temperature fields into steady (U,T) and fluctuating (u’,T’) components, 
the equation for the total kinetic energy of turbulence can be written 
(Lumley and Panofsky [1964]) 



fD — ’ 4 -1 p 0 • T* k hll \ + ^ ^fTo^O* • 6 Alt e - 4 

2 ^ rv T. ; 

® © © 

where Einstein’s summation convention holds, jx is dynamic viscosity, and 
T 0 the temperature of the undisturbed atmosphere. (^) is commonly called 
the production term; it represents the turbulent energy produced by fluc- 
tuations detracting energy from the mean motion by instability mechanisms 
as work=stress x strain. ^B) represents work done against buoyancy in 
terms of heat flux, IL = p T’u^ ; an upward heat flux acts as a source, 
a downward flux a sink of kinetic energy. (c) is the rate of working 



against viscous forces and represents dissipation of energy to heat ( € ). 
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The term -/>u’w’ from A is the Reynolds stress ( r ); it commonly 
represents momentum flux from the air to the sea surface (downward flux 
for > 0). This transfer across the boundary acts as an additional sink 
term as seen from the viewpoint of atmospheric turbulence. Similarly, 
it is apparent that the reaction stresses caused by the generation of 
waves acts as a source term for the generation of turbulence. Thus gen- 
eration of waves and generation of turbulence are inseparable. 

Horizontal (U) and vertical (W) velocities, temperature (T) , wave 
height ), Reynolds stress (r = -^u’w 1 ), and buoyancy • (/^ = w’T’) were 

measured or computed and related in time and frequency domains to depict 
cat’s paw structure. All data were taken at a single horizontal location, 
and w f and T were measured at only one height in air and water; it was 
therefore impossible to determine a meaningful quantitative energy bal- 
ance, although the qualitative interplay of production and stability terns 
was apparent. 



A.l. Individual Variables . Figures 2 through 7 (Figs. 2 - 
7) depict the imprint of three gusts. Fig. 2, from run 2, shows the 
fluctuations in 'rj , U, T^z , and /S for 40 seconds (4, 10-second time per- 
iods [T.P.]) starting just in advance of the leading edge of a gust and 
ending aft of the trailing adge. Fig. 3 shows samples of the computer 
generated (SC-4020) plots from which Fig. 2 was drawn (W is also shown). 

(The arrows and slant numerals on the x plot of Fig. 2 indicate the 2- 
second and 1-second segments chosen for detailed analysis.) Fig. 4 and 

Fig. 5 are similar representations of a cat’s paw from run 4. ,U, W, r , 

2 

a switch or intermit ten cy function of t ^ 100 dynes/cm , T, T 
r ^ are shown. The three non-consecutive T.P. are from before the 
gust, near the leading edge, and about 5 seconds later, still within the 
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gust. Fig. 6 and Fig. 7 are similar diagrams from run 6B. In run 6B 
there are three consecutive 10-second T.P. starting just at the leading 
edge, and one T.P. from 10 seconds after the gust. These gusts were 
selected to illustrate the full range of velocity and bulk stability 
values encountered (Table 1) . 

The basic features in all three runs areas follows. The leading 
edge of a gust is marked by rapid increases in both mean and turbulent 
wind velocities. Within a very short time (< 1 second), the air temper- 
ature drops as the boundary layer is mixed and destabilized by the tur- 
bulent stresses. T continues to drop as the winds increase; it starts 
to increase after U has reached its maximum value. There is commonly 
a secondary increase in U near the trailing edge that causes an additional 
drop of T. In run 6B, this second T drop is more rapid than the first, 
since the boundary layer has not yet recovered its initial static sta- 
bility and the stress fluctuations are consequently more effective in 
mixing. After gust passage, T rapidly rises to its initial value as the 
boundary layer restabilizes. There is no clear-cut relationship be- 
tween atmospheric variables and the low frequency wave field or water 
temperature. 

Both and yS traces are highly intermittent. There are bursts of 
high frequency, large magnitude fluctuations, separated by quieter 
intervals. Further, there are well-defined sub-bursts and intervals 
within the larger bursts. The r fluctuations are a maximum during in- 
creasing wind, but are not restricted to these periods. They are strong- 
est and contain most high frequency components not at the gust leading 
edge, but rather in the central bursts, after the boundary layer has 

been mixed and near the point of U . Bursts in r are generally asso- 
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Figure 2A 
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ciated with similar but shorter and weaker events in ; such 
coupled events are followed by T decreases. Bursts in also occur 
independently or with relatively stronger magnitude than the *r fluctua- 
tions. These events are followed by a levelling off or rise in T. The 
stress dominated fluctuations destabilize and mix the boundary layer, and 
buoyancy forces act rapidly to restratify it. The buoyancy forces are 
relatively stronger during run 6B, with high bulk stability. Here the 
sub-bursts are clearly evident. Comparison of the three runs, and of 
the two T.P. from run 6B, indicates that there are both stronger and 
higher frequency * and p fluctuations during higher winds. In the stable 
run 6B, fluctuations die out before those in ^ ; in nearly neutral 
run 2, even small z: anomalies cause variations in /6 . water fluctu- 

ations are also seen to be intermittent. They appear to occur during 

c 

atmospheric bursts and under breaking waves. T water fluctuations are 
very small, and strongly coupled to the wave field (i.e., they are due 
to motion of the thermistor through the vertical stratification); a 
small mean T water decrease results from mixing during the passage of 
the gust. (The spikes on the T^ ater traces of Fig. 5 are due to 
least bit noise in the digitizer.) 

The overall pattern indicates that fluctuations in the leading edge 
of the gust cause initial destabilization of the stratified boundary 
layer. This is followed by intermittent large magnitude instability 
fluctuations. * Buoyancy forces start to restabilize the boundary layer 
during quiescent intervals, and fully restore stability within 10-15 
seconds after the gust has passed. 

2 

The stress mean value s,<£/f oStt , are on the order- of 10 1 dynes/cm 
for run 2, 5 x 10° for run 4, and 10° for run 6B. <^> lOSec is approximately 
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_10 1 cm- ° C/sec for run 2, and +10 ^ for runs 4 and 6B. These values 

are consistent with the differences in mean conditions and comparable 

with the commonly assumed values used in dynamic modeling of the exchange 

—3 

processes (z = 1 dyne/cm 2 , H = 3x10 gm-cal/cm-sec) . The 1 instantaneous f 

3 2 

(250 hz) values commonly approach or exceed ~z ^ 10 dynes/cm and ~ 

2 

3x10 cm- 0 C/sec in the bursts. The strong high frequency fluctuations 

thus have amplitudes typically three orders of magnitude greater than 

the mean, or departures from the mean of 10 to 50 standard deviations. 

Similar types of measurements taken 5.6 to 22.6 m over a wheat field 

during unstable stratification (20 hz sampling rate) (Haugen, et al, [1971]) 

2 2 2 

yielded burst magnitudes of z - 10 dynes/cm , ~ 3x10 cm °C/sec. Burst 

duration was roughly 10 seconds, and U> M ~ 600 cm/sec. This data also 
showed the presence of bursts within bursts. McBean (1970) also found 
similar values over land, with a strong effect of stability (momentum 
transfer efficiency decreased and heat transfer efficiency increased with 
decreasing stability). The major qualitative patterns of intermittent 
gusts are thus similar over land and water. Quantitatively, the time 
interval selected for averaging is deterministic in estimates of instan- 
taneous fluctuation magnitude (Section C) . A 15-minute averaging period 
for the cat’s paw data would have increased the above estimates by at 
least an order of magnitude. The difference in magnitude between cat 
and land results is an indication of the significaice of stability and 
the enhancement of frequency and magnitude of stress fluctuations caused 
by the participation of a fluctuating boundary capable of supporting 
well-defined high frequency waves. 

Fig. 8 shows some typical histograms of the individual variables 
from the 10 second T.P. of all three runs. Fluctuation magnitudes and 
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Figure 8. 
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variance are much larger during than before or after the gusts; the peak- 
edness of the t: , ft , and w histograms is due to the high frequency 
fluctuations and consequent large number of axis crossings. U and T f 
distributions are more unimodal after than during the gust. There are 
few other basic differences between the T.P. f s for any one variable when 
the relatively long 10-second intervals are used. 

Fig. 9 shows typical autospectra based on the 10 second T.P. (6 degrees 
of freedom; 90% confidence limits are .63 x 10 n ^ to 1.74 x 10 n for a spec- 
tral density estimate of 10 n ) . With this comparatively long averaging 
time there are only generalized differences in frequency content between 
gusts and mean conditions. With the exception of the wave gauge data, 
the spectra are fairly white, with (note run 6B data) a slightly more 
rapid drop-off with increasing frequency in the T.P. after the gust. As 
also shown in the time plots, the high frequencies are augmented during 
the bursts. The differences in magnitude between gust and no gust con- 
ditions are apparent; the peaks of (f) and (f) are 3 decades 

lower after the gust. This basic similarity of pattern, but difference 
in high frequency content and major difference in intensity reflects the 
patterns displayed in Fig. 7 where the T and /3 scales have been deter- 
mined by the computer on the basis of the 10-second extrema. Wave propa- 
gation and the presence of long period wave groups prohibit any deductions 
from magnitudes of the wave spectra. 

Fig. 10 displays the changes in the 0.5 second vertical profiles of 
U from run 6B. It also shows <1 U > 0.5 sec and a measure of the turbu- 
lent kinetic energy production term. Profile perturbations, the devel- 
opment of unstable shear layers, and turbulent production occur predom- 
inantly during bursts. As with the r and^ fluctuations, production has 
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a range of i 3 to 4 orders of magnitude and is highly intermittent. The 
comparatively long averaging times CO *5 sec) are incapable of showing 
the extremely rapid fluctuations within the bursts; instead, the profile 
changes are similar to those that would be produced by a set of eddies 
propagating past the sensor array. Again, the maximum turbulence pro- 
duction occurs after the leading edge in T.P. 2, i.e., after the atmos- 
pheric boundary layer has been destabilized. The profile remains nearly 
steady, and production nearly zero, after the gust has passed and thermal 
stratification has been reestablished. Both stability and intermit tency 
are of major importance in gust time history. 

A sample time history of three cat’s paws thus indicates that gusts 
are highly structured. Averaging techniques such as spectra and histo- 
grams tend to obscure the significance of the microstructure. and yS 

fluctuations in bursts are at least 3 orders of magnitude above the 
mean. The localized instabilities also contain a larger range of fre- 
quency components than is the case in quieter periods and intervals. 
Production of turbulence and mixing in the boundary layer are strongly 
associated with the intermittent bursts; r fluctuations are largest 
after the mixing and destabilization associated with the leading edge. 
Participation of a large range of scales of motion and of the fluctuating 
water surface enhances the burst magnitudes. The use of long averaging 
periods obscures the details of the transfer processes that are dependent 
upon these strong localized instabilities. 

A. 2. Relationships between Variables . Relationships were 
depicted by joint histograms (JH) and time plots in the time domain, and 
by coherence (Coh) plots in the frequency domain. This duality of re- 
presentation helped to lend confidence to conclusions based on a limited 



39 . 



amount of data. Coherence may be defined as 

0,^(0 ■ k" ;ft 

** -i. us*>*< «?*>> 



It is represented in terms of its magnitude and phase, 

CoV^ (0 -- e l0I<<(0 . 

It was computed using fast fourier transform techniques. In the conven- 
tions adopted, ©^>0 means ^ leads ^ . Averaging was over 25 adjacent 

frequency bands for 10 second T.P., 9 for the shorter intervals/ As shown 
by McBean (1970), the ’spectral correlation coefficient 1 (or Coh) is an 
actual measure of transfer efficiency as a function of frequency (or scale). 



a. 10 Second Intervals . Fig. 11 shows JH(U,T I ) 

(with 10 second mean and trend both removed from T) from run 6B. Asso- 
ciated with the decrease and increase of < T > 10 sec is the boundary layer 
destabilization and mixing by the leading edge, violent mixing during 
the major bursts, and subsequent restabilization of a nearly constant temp- 
erature after gust passage. The top of Fig. 12 shows a similar pattern 
in the frequency domain using Coh r '^ from run 4. Before the gust the 
waves and temperature fluctuations (above 5 hz) are quite coherent and 
nearly out of phse. This is the pattern to be expected if the catamaran 
rides up and down through a thermally stratified layer that follows the 
pattern of the waves. (Ramzy and Young [1968]). Near the leading edge 

is slightly lower and extends to lower frequencies;©^ begins to 
shift, particularly in higher frequencies, as the atmosphere is mixed. 

Later in the gust the temperature field has been effectively inverted, and 
T and rj are more nearly in phase. The lower half of Fig. 12 shows Coh r ^ 
from the same T.P. Here is much lower but shows effectively 
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the same temperature inversion. 

Fig. 13 and Fig. 14 show JH(r,^ ) and Coh r ^ from run 6B; Fig. 15 
shows JH( r , ) from run 2. These diagrams show that the mixing and 

tumbling is effected by T fluctuations destabilizing the boundary 
layer, and that the following restabilization is due to the predominance 
of buoyancy forces. The change of the patterns with time is most clear 
in the more stable run 6B (note the change in scales); in the nearly 
neutral run 2 the fluctuations show less correlation, are more randomly 
scattered. Here the stress fluctuations predominate and remain intermit- 
tently large even after the gust has passed, while in run 6B the buoyancy 
forces rapidly debilitate the momentum flux as the gust departs. The 
coherence diagrams also indicate the phase reversals of stress and buoyan- 
cy fluctuations as the boundary layer is turned and mixed. The rapid, nearly 
linear shift of ©r/$ with f in the 9.5-17 hz range of T.P. 4 implies that 
fluctuations of similar frequency are ’slipping 1 past each other, or 
effectively that in the presence of strong stratification, weak z: fluc- 
tuations are unable to overturn and mix the stable boundary layer. 

Boundary layer stability is effective in retarding momentum flux 
by turbulent velocity fluctuations. This effect is particularly large 
near the leading edge of cat’s paws. Wave generation during high fre- 
quency instabilities will thus be less efficient during stable than 
neutral stratification, and delayed or reduced by the conversion from 
kinetic to potential energy involved in disruption of the stability 
structure. 

Fig. 16 and Fig. 17 show the JH and Coh represetations of the re- 
lationship between U at two heights during run 6B. As with the other 
variables, fluctuations are largest at and just after the leading edge, 
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with a collapse of the histogram after the gust has passed. is sig- 

nificant only at the DC level (T.P. 1-3) and the reversals of ®uo(P*C.) 
demonstrate the large variations in mean wind velocity shear that occur 
during the gust. The generally low level of ^ u u and lack of JH patterns 
over even this short distance are indicative of strong vertical struc- 
ture and the anisotropic character of the high frequency intermittent 
fluctuations . 

Fig. 18 and Fig. 19 are JH( z ) from run 6B and Coh from 

during and after the gusts of runs 2 and 6B and from all three T.P. of 
run 4. It is apparent that there is no effective evidence of strong 
wind-wave coupling during the gust when viewed from this averaging period. 
The effective generation mechanisms are obscured by the ’detritus f or 
by averaging over events of radically differnt character. If anything, 
z and y are more coherent in the 7-12 hz range (of runs 2 and 6B) 
after than during the gust; and in this frequency band the continuous 
linear ©ry drift is evidence that the fluctuations of similar frequen- 
cy slip by one another rather than strongly coupling. The turbulent 
stress-wave relationship that will be demonstrated has again been obscur- 
ed by considering averaging intervals much longer than the effective inter- 
mittent events. 

b. 2 and 1 Second Intervals . Two second segments 
containing primarily ’burst 1 or primarily ’interval’ type fluctuations 
were selected from runs 2 and 6B. The location of these 2-second T.P. 
is shown in Fig. 2 and Fig. 6 faith the exception of T.P. 1 from run 6B 
which was taken about 3 seconds prior to arrival of the leading edge.) 

A portion of the results is shown in Fig. 20, JH( Z , ^ ) from run 2, 
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Figure 19. 
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and in Fig. 21, Coh r y from run 6B. These figures show some overall 
similarity to those of the 10 second T.P.; some of the segments still 
contain both bursts and intervals, and this tends to partially blur the 
effects of the very high frequency fluctuations. 

Fig. 20 shows that the mixing and destabilization process is indeed 
associated with the bursts. Relatively weak, positively correlated fluc- 
tuations occur near the leading edge; mixing occurs during the next 
strong burst and is followed directly by negative correlation of -c 
and and consequent further destabilization of the boundary. With the 

boundary layer inverted the buoyancy forces work toward restabilization 
and rapidly degrade the strength of the momentum fluxes. Both the mixing 
and the restabilization occur rapidly; the intial boundary layer struc- 
ture is wiped out by the leading edge but is never firmly reestablished 

, . *y 

during gust passage. 

The wind-wave coupling begins to be apparent on this time scale. 

Fig. 21 shows that stress fluctuations are intially (before the gust) out 
of phase with the waves (7-15 hz range). They are almost totally inco- 
herent during the leading edge, i.e., during destabilization, then 

increases slightly in T.P. 3. This T.P. contains several well-defined sub- 
bursts and intervals so that the picture is confused, but is relative- 

ly constant at the peaks at 5.9 and 12 hz. During the major burst, 
the high frequency instabilities and high frequency waves serve to strong- 
ly and immediately couple the lower frequency fluctuations to the under- 
lying waves. The high frequency coupling is still not evident, but is 

high at and below the main wave period, with or the maximum 

stress occurring aft of the crests. In the immediately following inter- 
val at the D.C. end drops off slightly and the phase shifts; in 



52 




53 . 




54 . 



the higher frequencies is slightly higher but the propogating cap- 

illary waves now ? slip T with respect to stress fluctuations of the same 
frequency* 

In order to further decrease the blurring of burst and interval, 
segments approximately 1 second long were selected from within the 2- 
second T.P. (see Fig. 6). The results from run 6B are shown in Fig. 22- 
Fig. 24. T.P. 1 is prior to gust arrival; T.P. 2 just behind the leading 
edge; T.P. 3 during a short but intense fluctuation 2 seconds later; T.P. 

4 during the major burst and T.P. 5 from the quiescent interval 3-4 
seconds after T.P. 4. The pattern of events during this sequence may 
be described as follows. 

Prior to the arrival of the gust, t and y8 are nearly totally co- 
herent and in phase. Z fluctuations are small, loosely coupled to the 
main wave field, and out of phase with the high frequency waves. With 
the arrival of the leading edge of the cat T s paw the stress fluctuations 
increase in frequency and magnitude and the boundary layer is tumbled; 
the T ^ ft fluctuations decrease in coherence and* become nearly out of 
phase; little effect is noted immediately on the water surface, as the 
turbulent kinetic energy is utilized to stir the boundary; waves and 
are incoherent. After the intial burst the air is mixed and high fre- 
quency waves are generated, associated with the short, strong bursts of 
momentum flux. Z and j3> are again coherent, but out of phase, and the 
wind-wave coupling is apparently random. With the boundary layer de- 
stabilized, the magnitude of the stress fluctuations increases (T.P. 4) 
and the high frequency instabilities become strongly associated with the 
generation of surface waves. is quite low at frequencies above 15 hz 

with 6>r^ more variable than before but centered around O . 
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Increases considerably in the k - 2.8 hz range, but remains fairly 
low for lower frequencies. slips across part of the coherent ranges 

but shows a major reversal of trend in the center of the coherent band. 
During the quiescent period immediately after this event r fluctuations 
are much smaller and less strongly coupled to the short wave field, with 
slightly increased in lower frequencies, r 5 /6> fluctuations are 

again strongly coupled, and out of phase due to the mixing of the preceed- 
ing burst. This process is repeated during subsequent bursts, until 
the cat’s paw has passed the site and stability is reestablished. 

The change in character of the high frequency stress and wave pat- 
terns (Fig. 23) is strong evidence of the large growth rates associated 
with intermittent, large magnitude, high frequency instabilities. The 

generation of high frequency waves, as in T.P. 3 and especially in T.P. 4, 
* 

is indeed strongly coupled to the generation of turbulence. The large 
value of for lower frequencies in 2-second T.P. 4 indicates that 

the high frequency waves can rapidly and effectively couple the lower 
frequency waves and instabilities, and aid their -growth. It is obvious 
that little statistical confidence can be placed in low frequency measures 
based on such short time periods; nevertheless the change between time 
periods and averaging lengths, the similarity of this pattern in other 
bursts and gusts, and the enhanced high frequency in 1-second T.P. 4 

lend strong subjective credence to the hypothesis. 

Figs. 23 and 24 also clearly show how wave generation is delayed by 
the stability structure. Generation processes act only after the boundary 
layer has been initially destabilized. The continuous Gzrj slipping in 
T.P. 5 (across the full width of the coherent band) appears to be charac- 
teristic of ’detritus 1 , in this case, of the propagation and advection of 
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the waves and turbulence generated during the high frequency, localized 
instabilities. The slipping pattern is similar in T.P. 4 except for the 
10 - 18 hz range, where @ z rj levels off and reverses direction near the 
zero level. Most of the generation and momentum transfer occurs in this 
frequency range (see Fig. 27). 

Generation appears to be a very strong but selective and short- 
lived phenomenon. The time scales associated with the sub-bursts and 
instability associated transfers are very short, typically on the order 
of 1 to 1.5 seconds, which is also the typical period of the immediately 
underlying major wave field. The strongest coupling and high growth rates 
associated with the localized instabilities occur in a narrow high fre- 
quency band. Averaging times significantly longer than these scales can 
depict the general changes in atmospheric structure but fail to show 
the detailed nature of the coupling and transfer processes. The impor- 
tant scales are thus those of the waves, the bursts, and the internal 
fluctuations within the bursts. Momentum is transferred to the water 
at scales associated with these features. 

B. Interaction Between Scales . The microscales in both wind 

and wave fields account for only an insignificant portion of the total 
energy, yet as shown above they are the active participants in the local- 
ized high growth rate instabilities. They are effective in aiding the 
growth rate of lower frequency components and in coupling wind and wave 
fields, as was heuristically explained by Mollo-Christensen (1970, 1971) 
and as demonstrated in the time plots and coherence diagrams of the pre- 
ceeding section. Other evidence of nonlinear coupling between scales 
may be shown by an examination of the probability distribution of turbulent 
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velocity fluctuations and of the shape of the high frequency wave spectrum. 



has been shown to consist of a wide range of sizes of bursts or eddies 
of strong, high frequency fluctuations, separated by relatively quiet in- 
tervals, The cat’s paw itself is merely a burst on a lar^r time scale. 

The absolute magnitude of a single fluctuation is thus dependent not only 
on the instantaneous departure from adjacent values, but also on the de- 
parture of the fluctuation and its neighbors from similarly sized parcels, 
and on the departure of the mean magnitude of those patches from segments 
on yet larger scales; and so on up the range of scales to an averaging 
period covering the largest perturbations involved. Thus each scale of 
motion involves fluctuations whose magnitude depends both upon the 
strength of internal instabilities, and upon the forcing of and interac- 
tions between the scale itself and larger scales or larger anomalies. 
Following Cramer (1955), let the ins tantanious magnitude of a random 
variate be ; and let the change of this variate caused by a for- 

cing ^ be proportional both to the forcing and to the magnitude 
which had in turn been produced by the previous impulses. 



B.l. Scale Interaction in Velocity Fluctuations . A gust 
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which in turn, in the limit of small change, -** ~ 







For large N, the central limit theorem states that 



3- is normally 
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distributed. Hence 5; has the lognormal distribution, 







') 



a X =°> , 



where P ^ (x) is the probability 3 = x, CS " 2 is the standard deviation 



of 3 * E(x) = ^ § > , and a is an appropriate normalizing or cut-off 

value of 

The lognormal distribution can thus be expected to occur when a 
large number of impulses, or a large range of scales, interact to con- 
tribute to the increase of a variable at a rate proportional to its in- 
stantaneous value. For example, Rao and Narayanan (1971) have shown 
that the intervals between bursts in a turbulent boundary layer were log- 
normally distributed. This indicates that the scales of the bursts and 
intervals are determined by coupling with the scales of the boundary and 



of the outer or main flow as suggested by Mollo-Christensen (1970). 

Within a cat’s paw or gust, one variable which should act in a 
lognormal fashion is the turbulent energy dissipation, <£ . Dissipation 

is primarily dependent upon the magnitude of the -velocity fluctuation 



action between different scale eddies. The dissipation is due to the 
working of molecular viscous stress against this coupled rate of strain 
or deformation of the fluid particles. Provided that the range of scales 
involved in the deformations is sufficiently large, say including inter- 
action between the localized high frequency instabilities, fluctuations 
on the scale of bursts and intervals, and similarly between groups of 
sub-bursts and bursts up to the scale of the full cat’s paw, dissipation 
should be approximately lognormally distributed. 

Another variable which should be lognormally distributed, and for 



rate, which in turn depends strongly on the amount of ’rubbing’ or inter- 
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much the same reasons, is the magnitude of the squared time derivative 
of horizontal velocity. Time derivatives, like space derivatives, are en- 
hanced by the interactions between a variety of scales of eddies or in- 
stabilities. The interaction between high frequency velocity fluctuations 
and bursts has already been demonstrated in the preceeding section. Define 
the normalized squared magnitude of the velocity derivative as 



Chen (1971) has analyzed the distribution of with data taken 

from FLIP during BOMEX. He pointed out that squaring the derivative 
emphasizes the large fluctuations and minimizes the magnitude of small 
ones. (Gurvich and Yaglom [1967] showed that a lognormal distribution 

with large variance corresponds to a high probablity of values close to 

$ 

0 or much larger than the rms , but small probablities of intermediate 
values.) The natural logarithm of such a variable will have 60 
as ^ — ** O due to the ° ’tail 1 imposed by the rapid 

axis crossing of . Taking the average of a small number of 



terms this corresponds to averaging the fluctuations over a very small 
but finite volume, or numerically to limiting the analysis to scales 
slightly larger than the cut-off imposed by the sampling rate. This aver- 
aging will increase the tendency toward lognormality , provided that 
the number of scales interacting is sufficiently large and that the averag- 
ing time is short enough to include the important structure. Chen’s 
data were from a roughly 4-minute period at mean heights of 3 and 31 meters 
above the sea surface. The resulting distribution was almost identically 
lognormal for averages over 8 data points . 



( 3oi/ 3t y 




) wipes out this tail. In physical 
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Similar results from runs 2 and 6B of this study are shown in Figs. 
25 and 26. The burst/no burst/both character of the T.P. is noted in 
the figures. Fig. 25, with the data in histogram form (run 2, 2-second 
T.P.), clearly shows how averaging (over 5 points) removes the ’tail 1 of 



from both runs, 2-and 1-second T.P., are plotted on a normal probability 
scale. Purely lognormal distributions with comparable means and vari- 
ances are indicated by straight lines. 

Since only a relatively small number of points was used, the central 
limit theorem can not be strictly applied and none of the plots are 
strictly lognormal, though most of them show relatively minor departures. 
As shown best by the 1-second periods from the low stability run 2, per- 
iods of intense activity show a deficit of intermediate values, an excess 



of large and small ones. The quiet periods show the opposite effect, 
particularly in T.P. 6, taken after the passage of the gust (past the 
end of Fig. 2). In this very quiet T.P. the fluctuations are more nearly 
normally than lognormally distributed. The 2-second T.P. show these 
effects less clearly due to averaging over both bursts and intervals, 
but the tendency toward normality is clear in the quiescent T.P. 5. The 
run 6B data (high bulk stability) shows similar results, with intermediate 
deficits during bursts, excesses during intervals. One major exception 
is T.P. 2, from the leading edge of the cat’s-paw. It was shown in the 
preceeding section that in this time period the major effect of the 
stress fluctuations was the local tumbling and destabilization of the 
boundary layer. The excess of intermediate departures in this distribu- 
tion is further evidence of the initial effectiveness of the buoyancy 
forces in retarding momentum transfer and limiting fluctuation strength. 



the distribution. In Fig. 26, the cumulative probabilities of 
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The distribution of the squared velocity derivative is thus an ef- 
fective indicator of the range of scales involved in boundary layer tur- 
bulence. Interactions are strongest and involve the largest number of 
scales during the bursts. During quiescent intervals and particularly 
between gusts, fewer scales interact; coupling is weak or non-existant , 
and tends to have a normal distribution, indicative of a lack of 

interaction. 

It should be noted that for isotropic, homogeneous turbulence, i 
is related directly to dissipation by use of Taylor's frozen turbulence 
hypothesis, iL r - D i • Chen in fact regarded his data in 

these terms. Ruggles (1969) has demonstrated that the frozen turbulence 
hypothesis is quantitatively inapplicable in a case even less restrictive 
than this study; and indeed the above results indicate that it is not 
necessary to so restrict the assumptions on the nature of the flow in or- 
der to account for the lognormal behavior of fluctuation magnitude. The 
well-known intermittency of the small scale structure of turbulence, the 
consequent large dispersion of £ , and the dependence of the probability 

structure of € upon larger scales (Obukhov [1962], Landau [1959]) 
have been theoretically applied to the statistics of turbulence only for 
the locally isotropic, homogeneous case (Kolmogorov [1962]). The present 
development has avoided these restrictions and emphasized strong non- 
linear coupling between scales as the causative phenomenon. Further, 
it has been demonstrated that local departures from lognormality are 
also indicative of the relative strength of coupling of the turbulent 
components, and that many features of the boundary layer and transfer 
processes, such as the spacing of bursts and intervals and the frequency 
and strength of local velocity fluctuations, are controlled by this 
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coupling. 



B.2. An Example of Weak Resonant Wave Interactions , 



Nonlinear wave-wave interactions have been theoretically studied in con- 
siderable detail during the last 10 years (Kinsman [1965] provides a 
good summary of the major concepts). Energy exchange among gravity 
waves of different wave numbers is only seen when the perturbation ana- 
lysis extends to third order; it involves triads of primary waves inter- 
acting to feed energy into a fouth component. For waves in the cap- 
illary-gravity range, the analysis shows that interaction occurs at 
second order (as well as at third order); two waves interact to feed 
energy into a third (McGoldrick [1965], Simmons [1969]). This weak 
resonant interaction has been demonstrated in the laboratory (McGoldrick 
[1970]). * 

McGoldrick (1965) has discussed the implications of this energy 
transfer with respect to wind-wave generation. Of particular importance 
is the fact that the interaction times involved in the growth of the re- 
sonant wave are only 2-3 periods of the initial waves creating it. Be- 
cause of this short interaction time, the amount of energy that the wave 
field can receive from the wind at a given frequency is determined not 
just by breaking but also by nonlinear transfer to other components. 

Of particular interest is the 'self-resonant 1 wave, of wave number 



known as Wilton's wave; he noted the presence of a singularity at this 
wave number in his (1915) perturbation analysis. The general require- 
ments for resonance at 2nd order are k^+ ^2 - ^3 = ^ ~ ^3 = 




( X = surface tension/density). This is commonly 



and w; > J +T! W* ^ ^ 



where vj 



= 2.ivf AaJ k = w/h r lM 
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fulfills this relationship for k^ = k^ = ”1/2 k ^ > or ~ 



DC 



Simmons noted that waves shorter than are hydrodynamically unstable 
at 2nd order; further, k^ is the longest wave which can interact with 
two longer waves. Waves at k* are extremely unstable and can occur only 
as a transient state; they can not continuously exist as free surface 
waves and once generated will rapidly transfer their energy to the k^ 

C 2 



1 — 1 . 3 2 

waves. For clean water, JL » 75 cm /sec and k^ and k 9 roughly 



correspond to frequencies of 10 and 20 hz, respectively. 

Fig. 27 shows the wave autospectra from the 2-second segment ana- 
lysis of run 6B (80% confidence limits are roughly 1 decade wide). All 
5 T.P. display the peak of the major wave field at around 1 hz. T.P. 1, 
from before the gust, is rather smooth in the band above 3 hz. T.P. 2, 
near the leading, edge, shows a slight deficit around 9.5 hz, followed by 
a levelling off to 14 hz, a rapid drop to 19 hz, and a slight increase 
around 20 hz. During the ’active 1 T.P. 3 and the major burst in T.P. 4, 
spectral slope is less steep and the spectrum shows a slight bulge near 
10 hz (i.e., near the frequency where ©r rj leveled off and reversed 
slope); in T.P. 4, there is also a local maximum just above 20 hz. In 
T.P. 5, immediately after the burst, there is a significant deficit 
near 10 hz, and again a local minimum near 19 and. maximum just above 
20 hz. This sequence of events, particularly the change between T.P. 4 
and 5, is strongly suggestive of the picture given by McGoldrick. Dur- 
ing generation the k^ waves are not only present but near saturation 
because of the energy input from the wind; as soon as the forcing dis- 
appears, so do the extremely unstable k^ waves. Their energy is ’drained 1 
to higher frequencies. Fig. 28 shows a similar example from run 2, 
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where the mean wave field was significantly larger and the results con- 
sequently less clear. Here the major burst was in T.P. 3, with lesser 
activity in T.P. 4 and quiescent in T.P. 5.’ In T.P, 5 the spectral 
deficits occur at 10 and 18 hz, with the local maximum around 19.5 hz. 

The slight frequency shifts of the postulated and waves are 
readily attributable to contamination of the Charles River. The very 
fact that the predicted effects of weak resonant interactions stand out 
at all in other than a ’carefully tuned 1 laboratory experiment is 
strong evidence of their importance in wave spectral development. Ev- 
idence that perturbation analysis is suitable for predicting continuum 
nonlinear effects in major wave fields was recently provided by the 
bispectral analyses of Garrett (1970). Nonlinear transfer of energy 
is thus apparently of importance at all stages of spectral development. 

C. A Further Note on Time Averaging . The results of the preceed- 
ing sections have emphasized the importance of high frequency fluctuations, 
the intermittent nature of generation, and the coupling between scales. 
Averages over time scales comparable to the entire gust or to a combina- 
tion of burst and interval have been shown to obscure the selective 
high growth rates associated with the local instabilities, the effec- 
tiveness of wavelets in coupling the wind and the waves, and the nature 
of interactions between scales. One further example will serve to il- 
lustrate the critical importance of proper averaging in* the analysis of 
wind-wave generation data. 

Consider a typical cat’s paw from run 2. Typical mean velocity 
values for different scales can be very roughly taken as 
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Next, assuming rms turbulence levels of 10% (a very low estimate for 
a gust or burst), 3.75° C < 4.25°, 18 m/sec < u 'burst* ^ m /sec. 

As demonstrated in Section A, the important fluctuations in heat 
and momentum transfer and wave and turbulence generation are = w'T 1 

and T - —p u T w 1 (assuming symmetry in v f and V = 0) . Then 

U = < U > + u T , W = ^ W > + w 1 , T = T > + T f 

and since 

C W > and^wSare taken = 0 over whatever averaging period is chosen, 

u'w' = w ! (U -^U> ) ' and w ' T 1 = w 1 (T - C T > ) . 

Because of the impressed zero mean of W, the long term mean fluxes 
will be independent of the choice of U > and <^T> ; however, both the 
magnitude and the sense (sign) of the instantaneous fluxes (or of very 
short averages) are strongly affected by this choice. For the example 
above, choosing U = <U> g USt > T = <T>g^ s ^, and plotting a joint histo- 
gram of the resultant values of ZT and for the burst, the result 

will appear as in Fig. 29B (Fig. 29 is from run 2, 2-second T.P. 2); while 
if the means are chosen as the burst means, results are as shown in Fig. 29A. 
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When viewed against the general gust background it appears as if the ^ 
and fluctuations are negatively correlated (of opposing effect in 

turbulence production) and very strong; when viewed against the 2-second 
background, mostly burst, the f luctutations are relatively much smaller 
and show a weak positive correlation. Similarly, when the results are 
applied to the turbulent kinetic energy balance, the effects on production, 
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and hence on the stability structure of the boundary layer, are opposite. 



This is an extrememly straightforward but seldom applied result. The 
proper choice of averaging period must be made on the basis of the im- 
portant scales in the process under consideration. For wave generation 



and air-sea momentum transfer, these have been shown to be the scales of 
and within the bursts and the main wave field. Consideration of the wave 
generation process in terms of long term mean effects can totally obscure 
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the physical realities. 



Summary of Results and Conclusions 

The primary facts to emerge from this study may be summarized as fol- 
lows , 

1. Cat's paws are highly intermittent phenomena that them- 
selves are composed of localized bursts separated by quiescent intervals. 

The bursts contain high frequency fluctuations of stress and buoyancy 
3-4 orders of magnitude larger than that observed in the mean. The var- 
iability in time was shown in Fig. 2 - Fig. 7; representative magnitudes 
are shown there and in Fig. 8. 

2. The bursts are associated with the generation of high 
frequency waves and turbulence. This is best demonstrated by the se- 
quence of events shown in Fig. 23. Participation of the fluctuating 
boundary enhances the magnitude of the fluctuations , as shown by a com- 
parison of the results of this study with measurements over land. 

3. The high frequency instabilities not only have large growth 
rates on their internal scales, but also are associated with the growth of 
lower frequncy components. Both high and low frequency fluctuations are 
strongly enhanced in the bursts, as shown by the changes in magnitude and 
shape of autospectra (Fig. 9) and on the detailed SC-4020 plots. This is 
not only true for couplings within but also for couplings between media as 
demonstrated by the differences in coherence between stress and surface 
waves in the 2- and 1-second time periods (especially T.P. 4, Figs. 21 and 
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24). The bursts are also associated with the development of* shear layers 
in the wind profile and with large values of turbulence production (see 
Fig. 10), a further indication of the importance of localized instabilities. 

4. The boundary layer is mixed and tumbled by the instabil- 
ity fluctuations in the bursts. The effect of strong temperature strati- 
fication in the boundary layer is to delay the onset of .wave generation by 
requiring intial utilization of available turbulent kinetic energy for 
destabilization. Buoyancy forces are very effective in degrading the 
magnitude and effectiveness of momentum flux, and quickly cause restabil- 
ization after passage of the gust. 

The general mixing and boundary layer inversion patterns are shown 
in the joint histograms of U - T’ (Fig. 11) and (Figs. 13, 15) and 

in the phase reversals in the coherence of T f - and "c -/s (Figs. 12, 14). 
The mixing associated with the bursts, and the partial restabilization 
during intervals, is shown by the changes in correlation between ~c and 
/£ in Fig. 20. The delay of generation attributable to the utilization 
of kinetic energy for mixing is demonstrated in Figs. 22-24, where the 
burst at the leading edge (T.P. 2), uncoupled to the waves, has little 
effect upon the water surface but effectively inverts the boundary layer. 

The effects of stability have also been demonstrated by the ex- 
cess of intermediate sized velocity fluctuations in the leading edge 
(Fig. 26, runs T.P. 2), by the rapid rise of T after gust passage (Figs. 

11, 15), and by the temperature increases and decreases associated with 
Z and fluctuations of different relative strength (e.g., Fig. 6). 

5. The coupling between scales is strongly nonlinear, as 
illustrated by the lognormality of the distribution of squared velocity 
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fluctuations during bursts (Figs. 25, 26), the strong association of high 
frequency instabilities with large magnitude velocity, stress, and buoy- 
ancy perturbations (Figs. 2-7 and 9), and the association of burst length 
with the period of the major underlying wave components (Figs. 2, 4, 6). 

6. Weakly resonant nonlinear interactions are effective in 
transferring energy among capillary-gravity waves. This was demonstrated 
by the change in the and k^ components (Figs. 27, 28) between' ’genera- 
tion’ (burst) and ’detritus’ conditions. 

7. The important momentum transfer processes are of short 

duration, highly nonlinear, frequency selective, and strongly associated 
with intermittent high frequency instabilities and coupling between scales. 
This was demonstrated by the variations in coherence patterns between 

different averaging times and burst /no burst conditions (Figs. 21, 24), by 
the simultaneous occurrance of strong Reynolds stress and capillary wave 
activity , (Fig. 23), and by the selectivity of coherence and spectral change 
(Figs. 24, 27, 28). 

8. The differences between ’generation’ and ’detritus’ can be 

described by the changes in coherence between variables and in the proba- 
bility distribution of velocity fluctuations. Detritus coherence is char- 
acterized by low magnitude Oj ) and/or by continuous phase (© ) slippage 
in any coherent bands; during strong coupling and generation, T is high 
and & is nearly constant in the frequency bands that participate most 
strongly in the transfers (Figs. 19, 21, 24). Velocity fluctuations 
(measured as ^ ) in detritus have excesses of intermediate values and 

approach normal distributions; those from generation periods are more 
nearly lognormally distributed with excesses of very large values and 
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values near zero. These differences are due to the enhanced range of scales 
and strength of the couplings involved in wave generation and turbulent 
interactions. 

2 

The intermittency switch function (I = 1> ^ 100 dynes/cm ,1=0, 

2 

f < 100 dynes/cm ), shown only in Fig. 5 but computed in all runs and 
used for selection of the 2 and 1 second time periods, appears to be a 
suitable first approximation for a conditional sampling switch. When 
smoothed over a time period ^ 0.5 seconds it provided a reliable indicator 
of the presence of a burst and instability-associated generation. 

9. There is a rough, qualitative similarity between scales 
of motion (e.g., gust and burst), events of different magnitude (Fig. 7), 
and between turbulence patterns over land and water surfaces. The similar- 
ity lies in the presence of a large range of scales. The significant dif- 
ferences are caused by the strength of coupling between these scales; the 
larger the number of participating scales and media, the stronger the coupl- 
ing and the more efficient the associated heat and momentum fluxes. These 
features were demonstrated by a comparison of flux measurements over land 
and water, and illustrated in the time plots of the variables. The simul- 
taneous enhancement of frequency range and fluctuation magnitude in 
bursts is perhaps best displayed in Fig. 7, where the ’instantaneous 1 
magnitudes are based on the 10 second means. 

It was also demonstrated (Section C, Fig. 29) that the effect of ~c 
and fluctuations can be different when their magnitude and sense are 

evaluated over different scales; e.g., a burst may contribute to boundary 
layer instability as seen from the long term mean, while the local and 
immediate effect of the fluctuations is to temporarily reestablish a 
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portion of the initial temperature stratification. Averaging' over time 
scales other than those of importance in the momentum transfer process 
can lead to misinterpretation. 

10. Evaluation of generation in terms of the f mean f can be 
misleading and obscures the physics of the transfers. Since much wave and 
turbulence generation does not occur on the mean but in .localized, fre- 
quency selective, high growth rate instabilities Cas depicted in Figs. 23, 
24, and 27), and since the mean itself is a seldom if ever assumed state 
(as shown in the time series diagrams [Figs. 2-7] and the histograms 
[Fig. 8]), representation of the generation mechanisms only in these terms 
is unsatisfactory. 

On the basis of these results it is concluded that microscales play 
a predominant role in wave generation. High frequency waves and turbulence 
are intermittently generated during high growth rate, local instabilities, 
and the wavelets subsequently act as drag elements to couple the longer waves 
to the wind. Further, the high frequency instabilities aid the growth of 
similar, lower frequency fluctuations. This combination of events leads 
to spectral growth rates much larger than those predicted on the basis of 
linear and weakly nonlinear theories. 
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III. THE ROLE OF MICROSTRUCTURES IN WIND-WAVE SPECTRAL DEVELOPMENT 



Abstract 



Measurements of the wind, wave and surface current fields were made 
from a 38.5 meter spar buoy at the mouth of a tidal channel. Statistical 
analysis of velocity profiles showed that both air and water boundary lay- 
ers were predominantly logarithmic (neutral stability) , .and confirmed 
Ruggles T (1970) finding of anomalous behavior of atmospheric profile para- 
meters at critical wind speeds. A lack of simple relationships among pro- 
file measures indicates that wind-wave coupling can be only very poorly 
characterized on the basis of means, that the time history of spectral 
evolution is crucial to an understanding of the physics governing the mo- 
mentum transfer, and that more than one wave generation regime exists. 

A two-day history of spectral development, which did not include the 
major profile anomalies, was analyzed to determine the nature and scales 
of momentum transfer processes. Wave generation was found to occur pre- 
dominantly in short, intermittent, burst-type events, and to involve 
strpng nonlinear coupling of wind and wave fields with associated rapid 
growth of the high frequency components. These important events involve 
coupling between wave groups in wind and surface wave fields and strong 
interactions between the media and between frequency bands in both fields. 
Tidal currents modulate the wave field and cause imbalance of wind-wave 
equilibrium conditions. 



Introduction 

Microscale studies in cat’s paws have demonstrated the importance of 
intermittent, strong interactions between scales of motion and between at- 
mospheric turbulence and the sea surface in the generation of high frequen- 
cy waves. Similar types of interactions are also significant in the de- 
velopment of the major portion of the wind-wave spectrum. 

g 

In the open ocean, typical Reynolds numbers are on the order of 10 ; 
local wave generation occurs under the influence of surface drift, tidal 
currents, propagating swell, and stability, precipitation, wind speed 
and direction shifts associated with large weather systems. The large 
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number of participating scales and phenomena can all interact, with both 
strong and weak coupling, to have a significant effect on spectral de- 
velopment. It will be shown that there are separate regimes of wave gen- 
eration, and that the time history of events is of crucial importance in 
development of the wave field. 

Most previous open ocean studies of spectral development have dealt 
with mean or long time average wave energy growth and have examined the 
influence of only atmospheric conditions upon the wave field. Airborne 
measurements (e.g., Barnett and Wilkerson, 1967; Schule et al, 1971) and 
other laboratory and oceanic studies which extended over a range of fetch 
(e.g., Snyder and Cox, 1966; Sutherland, 1967; Mitsuyasu, 1969) have 
dealt extensively with the major features of spectral development with 
little or no experimental reference to the causative momentum transfer 
mechanisms. From these types of studies spectral growth rates, the fre- 
quency shift of the spectral peak, the existence of an equilibrium range 
for high frequencies, and the overshoot /undershoot phenomena have become 
well known features. 

Studies from stationary buoys or platforms (e.g., Seesholtz, 1968; 
Ruggles, 1969) have evaluated growth rates in terms of wind velocity pro- 
file characteristics with minimal consideration of the effects of the 
underlying current and subsurface boundary layer conditions. Such studies 
have shown that the atmospheric boundary layer is essentially logarithmic 
(at least for near-neutral stability) and that spectral growth rates are 
substantially larger than those predicted by the inviscid laminar Miles 
theory (Miles, 1957). 

Investigations of the direct interplay between wind and waves have 
concentrated on the transfer mechanisms suggested by this theory. One of 
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its critical elements is that turbulent and wave induced Reynolds stresses 
do not interact, and that momentum transfer occurs by means of an insta- 
bility in the mean shear flow. The momentum is extracted from the mean 
flow in the critical layer and transmitted downwards by wave induced Rey- 
nolds stresses. 

The recent work of Kendall (1970) demonstrated that a fixed wave 
pattern caused significant modulation of the turbulent structure* but 
that the induced fluctuations did not have a constant phase relation to 
the waves. Davis (1970) and Stewart (1970) theoretically and experimen- 
tally investigated the vertical structure of Reynolds stress. and found 
poor quantitative agreement with the wave-turbulence Reynolds stress pro- 
file predicted by instability theories. Wave induced Reynolds stresses 
accounted for only about half of the momentum carried to the wave surface 
by the pressure component in phase with wave slope. 

Ruggles (1969), Davidson (1970), and Shemdin and Lai (1970) have all 
demonstrated the existence of wave induced perturbations in the wind 
field and have shown that these wave induced Reynolds stresses do interact 
with turbulent stresses. There is thus substantial evidence that inter- 
action between atmospheric turbulence and surface waves is one of the cru- 
cial mechanisms involved in the transfer of momentum to the waves. This 
problem is being studied by Saeger and Reynolds (1971); their initial re- 
sults indicate that the inviscid model, by neglecting the important struc- 
ture of the critical layer, gives very poor predictions of the Reynolds 
stress and hence of momentum flux to the waves .. 

The cat f s paw studies have demonstrated the importance of detailed, 
near-surface measurements in determining the interactions among forces. 
Dobson (1969,1971) recently made the first successful detailed field 
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measurements of atmospheric pressure on wind generated waves. His im- 
portant results showed that the observed phase shifts were an order of 
magnitude larger than predicted by Miles' theory, and that energy and 
momentum flux spectra were strongly peaked at- or above the peak of the 
locally generated waves. Further, the inviscid laminar Miles model 
underpredicted increases in wave energy by factors of 5-8. These facts 
suggested a strong interaction between wave induced and turbulent Rey- 
nolds stresses. Dobson also noted that energy flux and wave supported 
momentum flux were negligible above relatively low frequencies, contrary 
to an assumption that much of the momentum transferred goes into short 
waves (Wu, 1969; and the cat's paw and present studies) but not contra- 
dicting the suggestion that short waves aid the growth of lower frequency 
instabilities by nonlinear coupling between scales of motion. Near sur- 
face pressure measurements over waves made by Elliott (1970) confirmed 
Dobson’s pressure shift findings. Elliott noted that pressure lagged 
the wave by about 135° during generating, and by 180° during non-genera- 
ting conditions. There was no vertical phase shift at heights above the 
crests. Elliott's work also suggested strong interactions between wave 
induced components and normal boundary layer turbulence. 

The basic conclusion from detailed investigations of near-navif ace 
structure and air-sea coupling is that turbulence and waves interact 
strongly and intermittently. Only a relatively small portion of the 
wave growth is due to processes which act on the ’mean', and these pro- 
cesses are themselves strongly dependent upon the presence and location 
of high frequency components and associated strong nonlinear coupling. 

The structure of the underwater boundary layer is also of major 
importance in transfer processes. It has been repeatedly demonstrated 
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that the magnitudes as well as the mechanisms involved in transfer of mo- 
mentum, heat, and moisture across an interface are radically different for 
flat plates or solid surfaces and fluctuating or moving surfaces (e.g., 
Chambers et al, 1970; Mangaralla et al, 1970). Over water, the processes 
are different for the cases of wind waves alone and wind waves with the 
added presence of swell and surface drift. 

Few experimental studies have considered the effect of the water 
near-surface layer motions upon wind-wave coupling, although the weak in- 
teractions between waves and currents or very long swell are well known 
theoretically (e.g., Longuet-Higgins and Stewart, 1960, 1961). The 
available experimental evidence on the structure of and interactions in 
the water boundary layer is somewhat contradictory. Shonting (1967, 

1968) noted that u-w phase relations were close to those predicted by 
classical theory, with no phase shift with depth and only minor interac- 
tion between wind waves and swell. Shonting, Stewart and Grant (1962) 
and Yefimov and Khristoforov (1969) have all found that wave stress is 
strongly concentrated very near the surface, and that essentially all 
dissipation occurs above the trough. The transfer between high and 
low frequencies caused by dissipation (breaking) and viscous decay of 
short waves (Longuet-Higgins, 1969b; Hasselman, 1971) is also a surface 
phenomenon, and involves non-linear coupling between scales. Shemdin 
and Lai (1970) found that the wind-induced drift in a wind/water tunnel 
had a logarithmic profile, and that the surface stresses computed from 
the air and water profile parameters were roughly equivalent. Intro- 
duction of mechanical waves disrupted the wind driven current. The 
couplings between the adjacent boundary layers and the interface are thus 
as yet imperfectly described. 
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This experiment was designed to investigate the nature and time 
scales of air-sea coupling during wave generation. It was conducted 
during a variety of wind, weather and current conditions. The analysis 
was divided into three segments. First, a study was made of the air and 
water boundary layers in terms of a logarithmic profile model. This 
section extends and complements the investigations of Ruggles (1970) 
and Groscup (1971), and describes the general environmental background 
of the more detailed analyses. Second, the time history of spectral 
growth during a two day period of increasing winds was investigated. In 
this section the general nature of coupling between wind, wave and current 
fields is described. Finally, segments of the two day record were studied 
in detail to emphasize the importance of strong, intermittent generation 
of high frequency waves and resultant wave- turbulence couplings that were 
demonstrated in the cat f s paw study. Generation ’events’ are described 
and their influence on the main wave field is displayed. 

The major objectives of the study were to demonstrate the wide 
range of interactions involved in wave generation and spectral develop- 
ment, and to emphasize the importance of high frequency ’microscales’ 
and intermittency in the momentum transfer process. 



The Experiment 

Measurements of wind and water speed profiles and wave height were 
taken from a stable, 38.5m spar buoy (Mollo-Christensen and Dorman, 1971) 
(Fig. 1). The buoy was slack moored in 31.1m of water (MLLW) and faced 
into the tidal current (NE flood, SW ebb, 1.5m/sec maximum current). 

The experiment ran from July through October, 1970. The site was at 
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igure 1. M.I.T. Air-Sea Interaction Buoy 
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the mouth of a tidal channel (Vineyard Sound, Massachusetts); minimum 
fetch to the S and SW (direction of prevailing winds and sea breeze) was 
50km. 

Instruments included 5 cup anemometers at heights of 2.4m to 11.4m, 

3 ducted-propeller (5cm diameter) current meters at depths of -1.5 to 
-4.7m, and two capacitance type wave gauges with 4.6m sensor length. The 
primary wave gauge was mounted directly over the current meters, 2.5m 
cross-stream from the anemometers. Thermistors were attached to arms of 
the velocity transducers at +5 and -2.5m. Buoy tilt was sensed with a 
pendulum pot. 

Signals from the transducers were mixed, telemetered to a shore sup- 
port facility on Cuttyhunk Island, Massachusetts and recorded in analog 
mode. Time, voice log, and 25kc clock signals were transmitted and re- 
corded for use in synchronization. Velocity information was in pulse 
format, wave gauge signals a variable frequency, and thermistor data a 
variable D.C. voltage. Data processing was done at M.I.T. and used a 
combination of analog and digital methods. Digitization rate for the 
time history analysis was lOhz; velocities were pre-filtered at 2.5hz, 
wave gauge signals at 5hz, to avoid aliasing. 

A total of 146, 10 minute ’mean condition 1 runs were made. There 
were 5 series of runs, with data taken once each hour for periods ranging 
from 6 to 53 hours, 109 of the runs contained data suitable for analysis. 
Winds and waves were from the S-SW during almost all runs, with 10 min- 
ute mean wind speeds (at 10m) of 1.5 to 12.5m/sec. Because of the large 
fetch the atmospheric boundary layer was nearly neutrally stable during 
all runs (AT air-sea ^ 3°C and -.02 < Ri(bulk^ ^ .05). Similarly, the 
water surface layers were generally well mixed down to a depth of at 
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3east 8-10m. 

The analysis was conducted in three parts. First, a distribution 
analysis was made of the parameters obtained from a comparison of velocity 
data with a neutral stability logarithmic profile model. Of the 109, 

10 minute anemometer runs, 96 (88%) had correlations 0.94 with the 
model. Of the 51 minimeter profiles used (from the longest nearly con- 
tinuous segment, runs 109-159, also used in the time history analysis), 

48 or 94% were logarithmic to within the required contraint; 49 of the 
anemometer profiles from the same period were logarithmic. Much of the 
analysis was based on 30 second averaging times. For the anemometers 
there were 2467 samples (81% logarithmic; 1057 samples with 91% loga- 
rithmic for runs 109-159), for the minimeters 1060 samples (85% loga- 
rithmic). Ruggles ( 1969, 1970) has discussed the significance of the 
correlation coefficient; Groscup (1971) has described the increasing 
tendency to a logarithmic profile with increasing averaging time for the 
runs discussed herein. Only the runs satisfying the Z 7 ^ 0.94 criterion 
were used in the histogram analysis, since the parameters are computed 
from the assumed model. This procedure is really only satisfactory for 
parameterization and description of the general background conditions 
and generation regimes, since as shown in the cat's paw study strong 
intermittent events are associated with large profile deflections and the 
development of unstable shear layers. The approach to a logarithmic pro- 
file with increasing averaging time is merely another indication that long 
term averages can obscure the physics of the transfer processes. 

The longest set of runs was then analyzed using standard time series 
techniques. The data is presented in terms of spectral growth and the 
related changes with time of atmospheric and water variables, and the 
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coherence-phase relationship between horizontal wind speed and waves. 
Spectral estimates have 6 to 20 degrees of freedom (dof); the 80% and 



are shown on the spectral growth and spectral' density estimate diagrams. 
Garrett (1970) has recently demonstrated the independence and normal dis- 
tribution of the real and imaginary parts of the Fourier coefficients im- 
plicit in the use of this measure for mean wave fields. The confidence 
limits are applicable only to each individual realization of the spectrum, 
and do not provide estimates of variance of the wavefield .for times longer 
than the overall averaging period; e.g., the 30 second spectra are strong- 
ly affected by the presence or absence of swell and major wave groups. 
Coherence magnitude ( t ) was tested on the assumption that the variates 
were random with a bivariate gaussian distribution (Amos and Koopmans , 
1963). (As demonstrated in the cat’s paw study and further displayed 
herein, this * a poor assumption for strongly coupled variables; un- 



fortunately it is the only readily available method). The results may be 
generally stated (for 5 minute time segments) as, at the .05 level of 
significance, the probability that the hypothesis = 0 is incorrect 



probability has dropped to 70%). 

Finally, several 10 minute segments from the start, middle, and end 



90% confidence limits of the ^ distribution (at <p - 10 n , n integer) 




for a measured ^ £ 0.5 is 97%. (For the 30 second segments, this 



of the set of 51 runs were analyzed using 30 second averaging periods. 
The data is presented in a similar fashion, with the added utilization 
of analog strip chart displays of ’events’. 
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Results and Discussion 



A. Boundary Layer Profiles . With the assumption that the averaged 
stress ( T ) in the near surface region is constant with height, a charac- 
teristic velocity of the turbulent boundary layer may be defined, 

U * = \f^ r J~ u w ! . 



This is termed the friction velocity. Using similarity arguments, for a 
neutrally stable atmosphere the only characteristic length is height. 

The resulting profile equation is 



B(.) - 



where r* 



is Von Karman’s constant - 0.40 and the roughness length, is general- 
ly not added to the height in the numerator (see, e.g., Lumley and Panofsky, 



1964). 



There is ^general agreement that this simple model gives a good 

description of the mean velocity profile in conditions of neutral stability 

over both land and water, but there is still considerable debate on the 

meaning of and relationships between the profile parameters (U. and z ) 

« o 

and other variables computed from the model. From the development of the 
model, and z q are independent similarity parameters characterizing 
the boundary layer; however, Sadeh, et al, (1971) have demonstrated that 
once the roughness elements are sufficiently high compared to boundary 
layer thickness, U. and z become local variables. 

Measurements of the profile parameters over waves have been re- 
viewed and discussed in detail by Wu (1969a, 1969b, 1970) and Ruggles 
(1969, 1970). Wu states that the wind stress ( T = ^ 0* ) is supported by 
form drag provided by wavelets superimposed on the larger waves; air flow 
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separates from the wavelets whose phase speed is less than U^. z^ becomes 
a measure of the height of these roughness elements (roughness height 
= 30 z q ) and is related to friction velocity by Charnock ? s relationship, 
zj ( U * 1 / <3 ) = constant. This is manifested by an increase in the 

drag coefficient (Cd z = r / P^i - ( ^ / ) ) with U up to a 'sat- 

uration 1 value, U = 15 m/sec. Above 15 m/sec the wind speed is larger 
than the average wave phase speed and Cd and z Q become nearly constant. 
Effectively, the Charnock relation implies that a single parameter can 
be used to describe the flow and thus the momentum balance of waves and 
atmosphere. The development of the profile model on the other hand as- 
sumed that and Z q are statistically independent. 

Ruggles showed that over waves, the wind speed is more appropriate- 
ly represented as the sum of a mean profile component (supporting atmos- 
pheric dynamics), and a term representing the correlation between and 
z^ (from wind-wave interaction). The 'correlation' wind component repre- 
sents that port(ion of atmospheric momentum required to support air-sea 
interaction dynamics. It is roughly 10% to 20% of U tota ]_ above 4 m/sec, 
and negative below 2 m/sec (indicating a minimum wind speed for generation). 
Ruggles also found marked changes in z^ and Cd , with excessively high peaks 
and excessively large variances, at around = 2, 4, and 8 m/sec. These 

anomalies correspond to the division points of the Beaufort scale, i.e. , 
to marked changes in the physical appearance of the sea surface. Similar 
anomalies are shown by Wu, who attributes the changes to the differences 
of flow in aerodynamically smooth, transition, and rough regimes. Hicks 
and Dyer (1970) found a simpler relationship, with Cd = .0011 and aero- 
dynamically smooth flow below 6 m/sec, and Gd = .0015 with rough flow for 
= 10 m/sec. They did not report any anomalies. 
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Additional evidence of anomalous wind speeds was given in the work 
of Groscup (1971) and Mollo-Chris tensen and Groscup (1971) (based on the 
data of this study). They further demonstrated a bi- and tri-modal dis- 
tribution of in the vicinity of the = 8 m/sec discontinuity, and 

a bimodal distribution of U. at or near the z and Cd discontinuities at 

* o 

4 and 8 m/sec. The anomalies are indicative of enhanced air-sea inter- 
action at specific wind speeds and of changes in the nature of wind-wave 
coupling. The large variances and bimodal distributions are evidence of 
enhanced nonlinear coupling at critical wind speed/surface ‘condition values. 

Fig. 2 shows the distribution (smoothed histograms) of 30 second 
average wind speed and of atmospheric logarithmic profile parameters from 
this study, both for all runs and for the subset 109-139 (those used in the 
time history analysis). Fig. 3 through Fig. 6 show selected joint histo- 
grams of the variables. 

U 10 was found from the profile model, using the computed values of 

and z . The drag coefficient, Cd, was determined for both the computed 

3 

10 m and measured 11.4 m speeds. is a measure of turbulent kinetic 

energy production, 

-Gv -- 

21 K i 



and Ou is the variance of velocity in the profile, 0\, c ( 0(0- o') 

The histograms show that a large percentage of the runs had very 

small roughness length ( (O * - 7 .5 or z 0 <\.8 y\6 Z cm; Values of 

_3 

z^ <1 10 cm were lumped into the In [z q ] = -10 category so that sta- 

tistics were not computed). Stress values are in rough agreement with 
the long term averages found by other investigators; the mean is largest 
for the continuous strong winds of runs 109 - 159, and decreases with 
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increasing averaging time. The values of Cd are comparable for the 
measured and computed heights; Ruggles (1969) computed Cd^ for 5 anemometer 
heights and 10 m, and found similar results > with as expected a slight 
decrease with increasing z. The Cd values are concentrated in the normal 
range for runs 109 - 159, but when all runs are considered together, the 
30-second and 10-minute means are excessively large and there are numerous 
exceptionally high drag values. These major, average changes in the drag 
effect of waves on the wind suggest the existence of more than one wave 
generation regime. <5 0 , Cc i lo > and (and thus and r ) all have 

decreasing mean values with increasing averaging times, indicating that 
profile distortion and resulting turbulence and wave generation occur in- 
frequently; the short, strong, important anomalies are obscured by long 
term averaging. 

Fig. 3 shows JH (U 1Q , U*), JH(U 1q , ln(z o ) ), and JH(U 1() , Cd ±Q ) for the 

30-second average profiles. In this figure, anomalies of U. ln(z ), and 

Cd^Q all occur at roughly 2.5, 3.5, 6.4, 7.7, 10, 11, and 12.2 m/sec (they 

are particularly striking at 2.5, 3.5, and 7.7 m/sec, i.e., near those 

found by Ruggles). The anomalies are in the form of exceptionally large 

variances and wide departures of the interval means from any postulated 

simple fit of the data. The values are in general quite scattered and 

2 

poorly represented by means. Also shown is JlUU^ /g, K) where K = 30 
is height of the roughness elements (wavelets) . as postulated by Wu. It 
is obvious that a simple relation such as suggested by Charnok and found 
by Wu does not hold here. There are bimodal distributions and anomalies 
in both directions. A single parameter can not characterize the profile, 
because of the interplay between atmospheric dynamics and the wave field. 
None of the data can be fit with a linear relationship; in particular. 
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Ruggles's suggested = 0.04 U^q intersects the minima of the mean of 



U* on u 10 » 



Fig. 4 shows Cd.^) com P ute d using the formulae of Wu (1969b) 



for 1 breeze 1 and * light wind* regimes, 

z 

U £ ^ “Aex * 



Ci , 

to 



J..( ^.lu» OQ- ) 



A*/ . ) 

.OlSC. Ua /<} / 



,*/ <ru<r 15 M / 



The anomalies at 2.5 - 3.5, 7.5, and 12 m/sec again appear, but much less 
strongly than in direct computation from the profile model; in addition, 
the mean value of Cd computed using this technique is excessively small. 
Large scatter was found in joint histograms of other relationships sugges- 
ted by Wu. The failures of these postulated simple models on the relative- 
ly short (30 sec) averaging times indicates that the dynamics of wind flow 
over open ocean wave fields at large Reynolds numbers can be only poorly 
described in terms of long time averages, even when the significance of 
wavelets is considered. Wu ! s contention that the small surface components 
support the majority of wind stress is not disproven by these data; the 
results merely indicate that the dynamics of the wind-wave coupling is 
considerably more complex than that caused by just form drag and air flow 
separation. 

Fig. 4 also displays 30 second JH( r , ln(z Q ) ), JH(Cd^Q, ln(z Q ) ) and 

JH( /k? )• The T -ln(z^) relationship shows a distinct bimodality, 

with generally increasing stress with roughness length and stress maxima 

for z q ^ lew (ln(z Q ) = 0). High stress occurs predominantly over large 

’roughness elements’, but only under certain conditions as yet undetermined. 

2 

The Cd- A = (U./U ir .) vs ln(z ) plot is very closely fit by the profile 

ID * 10 o 

relationship, (t*) * 

(using K = 0,40; Businger, et al, 1971, have recently noted that = 0.35 
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is a better value of Von Karman 1 s constant for neutral stability over a 
land surface); this is -merely a check on the consistency of the logarith- 
mic profile assumption and on the continuity of data from different time 
periods, and serves to emphasize the significance of the anomalies of 
the other plots (i.e., they appear not to be due to instrumental or analysis 
peculiarities, and the 0.94 correlation is a sufficient criterion for 
assuming a logarithmic distribution). The cr x u - 0* 5, relationship 

shows that turbulent production increases with increasing fluctuations, but 
with a very strong bimodality, indicative of the presence of more that one 
production-generation regime . 

Fig. 5 shows JH(U^q, ln(z Q ) ) and JH(U^q, Cd.^) for averaging periods 
of 2 and 10 minutes. The basic patterns of anomalies at specific values 
of persist at the longer averaging times, but with less variance. (The 

10-minute averages have been compared by Groscup [1971] to those of Ruggles 
[1970] [on the basis of loglO (z q )]). The basic patterns remain, irrespec- 
tive of averaging time, but the decreasing intensity of the anomalies shows 
that, as in the cat’s paw studies, the strong and intermittent events are 
blurred or obscured by long time averaging. Even with these long averages 
and highly logarithmic profiles, it remains clear that the ’mean 1 is a 
very poor characterization of physical reality. The time history of events 
and the regimes and instantaneous characteristics of wind and wave fields, 
are highly significant in determining the nature of the boundary layer dy- 
namics . 

This point is further emphasized in Fig. 6, which shows JH of with 
U^, In (z q ) , and Cd^Q for runs 109 - 159, and JR(U^q, U^) for all runs for 
the two subclasses | (4. < ^ <T )o The JH of 

runs 109 - 159, characterized by strong and increasing wind and a large wave 
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field, occupy concentrated segments of the overall diagrams. Signifi- 
cantly, most of the anomalies are missing, and as previously noted Cd.^ 

-.3 

for these runs was 0.9 x 10 lower than the overall average; the value 
for these runs is in rough agreement with that’ found by most investiga- 
tors. This set of runs apparently satisfies the criteria for a single 
or homogeneous wave generation regime (the bimodalities of the other rela- 
tionships also disappear or decrease when these runs are plotted separate- 
ly) . The division by mean velocity derivative magnitude shows that the 
parameter relationships do vary with time history, and that (as in the 
cat's paw study) higher wind speeds are in general associated with stron- 
ger fluctuations. 

Fig. 7 shows examples of the results of the JH of the parameters of 
the logarithmic underwater profile (runs 109 - 159). 30-second mean 
current speed (at 10m) was approximately 10 cm/sec, with a range of -3 
to 27 cm/sec. The diagrams show quite consistent relationships among 
the profile values, with no indication of the anomalies or bimodalities 
found in the atmosphere. values are negative, as the wind drift is 

strongest near the surface. Joint histograms between air and water pro- 
file parameters from this set of runs failed to show any overall sta- 
tistical relationships; this is not surprising in view of the fact that 
most of the momentum from the wind goes directly into the wave field 
(Stewart, 1967; Dobson, 1971), so that the effect on the surface drift 
current is not immediate or else (see section B.2, Figs. 28 and 29) causes 
the short average profiles to become highly nonlogarithmic. 

In summary, the mean velocity profiles in both air and water are 
predominantly logarithmic during conditions of neutral stability. Ano- 
malies in the atmospheric profile parameters demonstrate the importance of 
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time history and the existence of more than one wind-wave generation 
regime, with enhanced activity at certain critical wind speeds under as 
yet unknown circumstances. Anomalies standout most strongly and the 
air boundary layer becomes much less logarithmic at short averaging 
times, indicating the importance of short, intermittent events in the 
development of unstable shear layers and associated instabilities. The 
profile can not be well represented by a single similarity parameter, and 
there is no simple relationship between wind velocity and the profile 
measures. Attempted use of such relationships for characterizing the 
wind regime and predicting wave growth can lead to significant errors. 
Consideration must be given to the time history of generation and to the 
importance of intermit tency, i.e., to the enhanced air-sea coupling and 
momentum transfer at certain wind speeds or in short events. 
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B. Time History 

The 51, 10 minute/hour runs (consecutive except for 1 hour gaps 
between runs 20 and 21 and 42 and 43) selected for detailed analysis were 
taken from 21 to 23 October, 1970. The wind direction had switched from 
N to S 30 hours prior to the first run, as a large high pressure system 
(> 1040mb at the center) passed the site moving NE. Winds increased and 
weather deteriorated throughout the three day period, then winds slacked 
and weather cleared on October 24. Analysis of current records and of the 
tide record at Wood’s Hole, Massachusetts, showed that the surface drift 
continued NE, into Vineyard Sound, after the NE flood tide on the morning 
of October 21. The buoy thus faced into the wind and waves throughout 
this set of runs. This was confirmed by visual observation. 

B. 1. 5 and 10 Minute Averages 

Figs. 8 and 9 depict the time changes of the 10 minute 
average horizontal wind speed at 10m (computed) and 2.4m (measured) and 
the wind extrema at 2.4m. The figures also show current speed at Z = 
-1.52m, and stress ( r zp U* ) computed from air and water profiles. Short 
lines connect the averages of 2, 5 minute periods from each run and give 
an indication of the variability on shorter averaging times (the large 
offsets in Runs 111 and 142 are due to noise; those in 120 and 157 are 
real; 2.4ra velocity extrema from runs 129-133 were unavailable). 

There are strong periodicities in all of the variables. Mean wind 
shows an overall increase from 8-12m/sec; winds increased slowly through 
run 132, then increased rapidly to lOra/sec, stayed near that level for 10 
hours, decreased and remained moderate for a similar period, and then 
again increased rapidly starting around run 152. Superimposed upon this 
long modulation is a shorter but strong periodicity of roughly 2-4 hours. 
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Velocity shear and extrema also show an overall increase with time, with 
superimposed shorter periodicities. The water current has a strong semi- 
diurnal tidal component with alternating times between maximum flood of 
5 and 7 hours; the current maxima around run 139 is strongly suppressed. 
Some of the wind fluctuations occur with the tidal frequency and its 2nd 
harmonic; local wind maxima frequently occur during maximum flood. 

The effect of the tidal modulation is vividly displayed in the 
water boundary layer profile as an excess during ebb, when the wind and 
tidal components are opposing and relative wind to water speeds a maximum. 
The short 2-4 hour fluctuations of the wave field are also apparent in the 
water stress. Air stress has a diurnal periodicity and a small 2-4 hour 
component as well as local maxima during the two major wind increases; it 
also temporarily increases in runs 141 - 146 when the winds decrease ra- 
pidly. Mean stresses are thus largest when there is a significant distur- 
bance of the equilibrium balance between current, wind and wave fields. 

The stress values from air and water are in rough agreement. Differences 
may be partly accounted for by the momentum carried by propagating swell 
and wind waves, and by the momentum and equilibrium inbalances caused 
by the tides. 

Mollo-Christensen (personal communication) has noted that- except for 
the tidal variations, 10 minute mean air and water velocities are compar- 
able when normalized by their dynamic viscosities, i.e., . 

This is shown' in Fig. 10 (from Mollo-Christensen and Groscup [1971]). 

This relationship can be accounted for by assuming that the main shear is 
in the viscous sublayer, so that 
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If •=- , e.g., if the sublayer thickness in both media is controlled 

by the same elements (say the wavelets), then since the stresses are in 
rough balance the boundary layers are shear coupled and related by the 
ratio of their dynamic viscosities, U a r Tp ^ • The long term ave- 

rage values thus again support the assumption that the major air sea coup- 
ling can be simply described, in spite of the known intermittent discon- 




Figure 10 

term variations of air and water stress. 

Fig. 11, Fig. 12 and Fig. 13 show the variability with time of the 
wave field parameters. Fig. 14 shows the corresponding wave spectra, and 
Fig. 15 the development of selected frequency components (cross sections 
of Fig. 14). 80% and 90% confidence limits at = 10^ are shown by ticks 

inside the axis of Fig. 15. E is mean energy density, E = P^^\ Z 
and momentum density M = E/c (c = wave celerity) computed by spectral 
summation. c ) and M (§ MA< ) are respectively the celerity of the 

wave of maximum spectral density and the momentum computed by assuming 
that the energy density is concentrated at that one frequency. The 
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information in these figures may be summarized as follows. 

During the two-day period the mean energy density of the wave field 

increased by nearly an order of magnitude, from 2 x 10“* to 17 x 10”* ergs/ 

2 

cm . The period of the wave of maximum spectral density (T) shifted to 
lower frequencies, from 4 seconds at run 109 to >■ 8 seconds at the end. 
These changes can be clearly seen in Fig. 11, Fig. 12 and Fig. 14. The 
changes in T and E are slow and comparatively steady for the first 23 
hours. E then increased rapidly in conjunction with the rapid wind in- 
crease that started around run 132. Fig. 14 shows that the spectral peak 
turned sharply toward lower frequencies during this velocity increase; 
after run 132 the separation of the spectrum into two peaks (swell and 
locally generated waves) becomes more distinct and the peaks continue 
to separate as the swell peak shifts to still lower frequencies. The 
equilibrium range above the local peaks has a very rough f ”* slope 
(Phillips [1958]) with occasional local maxima in the range of 0.5 to 1 
hz. The distinct separation between local and swell peaks was wiped 
out during the rapid wind increase of runs 148 and 149. 

All of the figures show that the growth of the wave field, in terms 
of both mean energy density and of individual spectral components, is not 
steady. There are large, significant variations of E, T, normalized par- 
ameters, and of (-f ), the spectral density at selected frequencies 

(Fig. 15). E (1 Hour) (Fig. 13) becomes as large as 4- ,1% of the hourly 
average E during later runs. The hourly changes, both growth and decay, 
are largest when the wave field is most energetic and during major wind 
speed changes. 

Fig. 15 demonstrates that even the higher frequency components (.6310 
and .3550 hz) , which are always above the spectral peak and in the f equili- 
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brium* or saturation range, have significant (at the 90% confidence level) 
increases and decreases of energy content on these one hour averaging 
periods. The .2515 hz component, which was at the spectral peak at the 
start of the runs, appears to ’overshoot 1 in runs 110-112. Its spectral 
density then falls off to a lower, ’equilibrium 1 value during the slow 
growth period that extends through run 132. After that run, the 4 second 
waves again show major growth, ’overshoot’ slightly at the same spectral 
density they had earlier reached, and then settle into a second, higher 
’equilibrium’ level. In runs 154 - 156 their energy content increases to 
still higher levels. The wave field can apparently support more than one 
overshoot and saturation range under different wind conditions. 

The .1583 and .1254 hz components both grow until around run 139, 
when the higher of the two frequencies reaches a quasi-steady state level. 

The most rapid growth of the low frequency waves occurs in the rapid wind 
increase starting at run 132. The major wind decrease of runs 142 - 143 
is not accompanied by significant decay; the <^(f) decrease, and E(1 hour) 
for this period, are no ]arger than those found in the minor fluctuations 
during growth periods. The fluctuations of the individual frequency 
components is not in phase; some portions of the spectrum grow, others 
decay, in what appears to be a random fashion except during the major 
wind increase at run 132, which is followed around run 135 by growth of 
all but the highest frequency components. The locally generated waves, 

.2515 and .3550 hz, show similar periodicities to those of the local wind 
field* The individual swell and high frequency components are essential- 
independent of these short term variations, while the total wave field - 
E, T and M - does show the effect of local generation. 
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A comparison of the plot of U^/C ) with those of E and M 

shows that spectral development is most rapid when the wind is moving fas- 
ter than the waves of the spectral peak. Some growth still occurs even 
when the waves move faster than the wind (where Miles theory is clearly 
inapplicable due to the absence of a ’critical layer 1 ). The general 
range of U /C ( ) is quite low for the entire set of runs. These 

features indicate that swell was propagating into the area, but do not ac- 
count for the lack of local decay of low frequency components between runs 
137 and 158 when U /C ) was continuously less than 1. It may be 

hypothesized that energy was being transferred to the low frequency waves 
by nonlinear resonant interactions within the wave field, (e.g., Garrett 
[1970]) or that wind generation occurs at these frequencies during brief 
but intense periods when U was larger than C ( ^W v ). 

The plot of C ( § x )/U^ shows a strong diurnal periodicity, and 

lags the similar variation in I .by roughly "H/2. There is no ap- 
art 

parent relationship between this parameter and U^/C ( max ) > or between 

it and energy density growth rates. 

The normalized wave momentum plot indicates that only 60-80% of the 

total wave momentum is carried by waves at the spectral peak; most of the 

remainder is carried by waves at higher frequencies (lower phase speeds). 

Fig. 16 shows the momentum spectrum (5 minute averages, 6 dof) from several 

runs (ordinates have been normalized to the spectral maximum). These 

plots show the overall increase in M and M , the shift of M to lower 

max max 

frequencies, the bi- and trimodality associated with the separation of 
swell and local wind wave peaks, and the occasional very strong concentra- 
tion of momentum in a single swell frequency. The momentum spectra have 
roughly the same shape as the power spectra, but tend to emphasize the 



116 . 




117 . 



concentration of the wave momentum in one or more narrow frequency bands. 

The 10% momentum levels — the range of f for M(f) .1M — were 

b max 

approximately .15 to .55 hz through run 131, and then decreased slowly to 
roughly .09 - .40 hz. Wave momentum was indistinguishable from noise 
level above 1.40 hz near the start of the runs, and this limit decreased, 
starting at run 132, to approximately 1.00 - 1.20 hz. Most of the mo- 
mentum spectra had several discrete peaks, and there was usually a consid- 
erable difference in the minor features of the 2,5 minute segments in 
each run. Current spectra at -1.5m (not shown) had maxima at both swell 
and locally generated wave frequencies, and were generally coherent with 
both of these wave field components. 

Coherence between horizontal wind speed at 2.4m and waves at the 



0.65 for 50%, and ^ 0.50 for 97% of the 102, 5 minute samples. 
Phase ( ) at that peak was IT for 82% of the coherent runs, and 

> ^ 2 radians for all of them. As noted by Seeshdtz (1968), wind 

speeds are maximum over the wave troughs. 



the peak to lower frequencies with time, and its bifurcation into local 
and swell segments, are clearly evident; there was a secondary maximum 



the samples after run 131. In the majority of these cases (note run 




measured on the 5 minute sub intervals) was 





137-2) there was a decrease in V between the two peaks, with a shift 
in from Tl to p/t. and then back to V . Winds were generally in- 

coherent with the waves ( To* statistically indistinguishable from 0) 



above the local wave peak on this averaging time; one notable exception 
is segment 157-1; this run contained a very strong generation event and 
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is discussed in detail in section B.2. In roughly 30% of the runs , $ Orj was 

fairly high at the D.C. level; in all of these cases, . With the 

exception of the occasional D.C. peak, Xov> was low for £ < f J 

( ^ max 

Wind spectra did not consistently display distinct peaks at the wave 

-c 4 

spectral maxima (Fig. 18, 80% and 90% confidence limits at = 10 and 
= 10 5 are shown on the "3? axis of run 142). They were fairly white 
( < 1 decade dropoff in the . 2 to 2 hz band). Nevertheless there was 
nearly always a spectral increase at or around f £ max > and an cases °f 
a double-peaked coherence the wind spectra showed a local minimum between 
the peaks. The examples show (f) and <3? ou (f) for the same run seg- 

ments as in Fig. 17 (they are rough, 6 dof, and the jaggedness is probab- 
ly not significant). 

From these facts and figures the following inferences may be drawn. 
Wave generation and spectral development is strongly non-s tat ionary . There 
are significant, rapid changes in spectral growth and these changes are 
poorly represented by long term means; no single parameter or simple com- 
bination of parameters can describe the process. Air and water boundary 
layers (both velocities and stress) are viscously shear coupled. Modu- 
lations at the tidal frequency and its harmonics ( or of the period 
ATide/ A t) appear in the wind and wave fields; changes in the tidal flow 
disrupt the equilibrium between wind and wave fields and cause modifica- 
tion of momentum advection patterns. 

Even though the high frequency portions of the wave field are near 
saturation and account for insignificant portions of the spectral density 
and momentum, they show significant , rapid (hourly) variations and can 
account for spectral growth when the relationship would normally 

predict major continuous damping. Measured growth rates (long term 
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averages) of these high frequency components are inherently incapable of 
displaying major changes, since on the mean these waves are near satura- 
tion, and also because the very high frequencies are continuously being 
generated by the larger waves (Longuet-Higgins , [1963] ; Crapper [1970]). 

The importance of these high frequency waves lies in their ability to 
couple the wind with larger waves and to assist the growth of low fre- 
quency instabilities by nonlinear coupling between scales. It will be 
shown in the next section that these components have extremely rapid 
growth rates in intermittent generation ’events’, but again these impor- 
tant physical momentum transfer features are obscured by the averaging 
process . 

The influence of waves on the turbulent atmospheric boundary layer 
is strong at both the swell and the locally generated spectral peaks. 
Furthermore, the 2-4 hour periodicity of wind and both the locally gen- 
erated and overall wave fields are strongly related. Changes in measures 
of the wind and wave fields (e.g., E (1 hour)) are largest when the mag- 
nitude of the measures is largest, similar in nature to the pattern of 
higher frequency components during strong gusts that was noted in the 
cat’s paw study. 

The major conclusion from this section is that study of the time his- 
tory of generation and spectral development is essential for an understan- 
ding of the physics involved. Wind, wave and current fields show evidence 
of strong interrelationships that are poorly represented by parameteriza- 
tion of the mean. The rapid variations in E, E (1 hour), U and £ all 
indicate that generation is predominantly an intermittent process, char- 
acterized by rapid transfer of momentum and followed by weaker couplings 
and nonlinear interchange between spectral components. 
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Figure 21. 
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B.2. 30 Second Averages . 

Several of the 10 minute runs were chosen for de~ 
tailed analysis to help determine the nature and scales of intermittency 
and the resultant momentum transfer. Analysis and data display were simi- 
lar to those used for the 10 minute runs. 

F-?g- IQ a nd F xg. 20 depict the changes with time of 30 s econd 



U (cm/sec) 
1200 - 



900 



600 - 



“80 




80 “80 0 80 
7 } (cm) 

Run 110-10 Run 159-3 



JOINT HISTOGRAMS 

U ' 7 ] 

Samples per Interval 

I “ 2 
- 3 “ 4 

• 5-6 

• > 7 



30 second samples 



Figure 22 

averages of variables in the air and water for runs 110, 127, 157 and 159. 

Fig. 21 shows E (30 sec) and C ( <R ) /U for the same runs. The major 

max * 

features of the 10 minute runs also stand out in these displays of short 
term change. There is a general increase in U, E, M, T, wave height and 
spectral density and momentum extrema. Stress computed from the air and 
water profiles are in rough agreement, although they show different perio- 
dicities. The short term changes are largest and most rapid when the mag- 
nitudes are largest; changes in U are paralleled by increases in E, though 
not as clearly as in the longer averages because of the influence of wave 
groups (these had a roughly 40 second periodicity for later runs, and also 
significantly affected the T w measurements). U is again normally over 
the troughs, but with much scatter; this is illustrated in Fig. 22, JH 
(U ? ., YJ) from 2 short segments, and is more clearly shown in the coherence 
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diagrams (Figs. 29 and 30). 

In these runs, the waves of the spectral peak carry roughly 80% of 
the total momentum density (again implying that a significant portion is 
carried by waves of higher frequency, or slower phase speed) , while the 
waves of the momentum peak (when different from the spectral peak) over- 
account for the total momentum. U^^/C($> ) is in the range of 1 - 2, 

10 max 

and is in fact largest during the earliest run, due to the absence of 
swell. This parameter shows no significant peaks at the points of wind 
and E maxima; this may be interpreted to mean that the wind is strongest 
over the largest wave components, indication of selective coupling. 

C ( )/U, shows no distinct relation to spectral growth or decay; 

E (30 sec) has comparable variations in runs 110 and 127, and the para- 
meter is similar in magnitude in runs 110 and 159 when the energy growth 
patterns are radically different. 

For these short averaging times E(30 sec) is 1% to 10% of the mean 
energy density for the corresponding 10 minute periods. The rapid changes 
in spectral shape and content are well shown in Fig. 23, from run 157. 
Cross-sections of this diagram are. shown in Fig. 24 (n.b. 80% and 90% 
confidence limits at 10^ on axis). It is immediately apparent that 

the lower frequencies (referring here both to swell [.1254 hz ] and to 
local wind generated waves [.2515 hz]) are relatively constant, with 
variations primarily attributable to propagating wave groups (this was 
checked by examination of analog strip chart displays of the filtered 
signals). The change in the high frequency components, however, is 
very dramatic, with 30-second variations of up to 2 orders of magnitude 
(3 orders of magnitude for the 4.46 hz waves. Fig. 23). The change of 
shape of the ’equilibrium range’ in Fig. 23 also demonstrates this 
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variation; there is a repeated pattern of large excesses in specific high 
frequencies followed by a filling out of the spectrum. The subjective im- 
pression of high frequency excesses climbing the slope of the spectrum 
and transferring their energy to lower frequency components is hard to avoid. 
Fig. 24 also shows that the excesses do not always occur simultaneously in 
all frequency bands and that on the 30 second averaging period the high fre- 
quency growth does not always correspond to the presence of low fv-equency 
peaks. Further, intermittent high growth rate events are of short duration, 
typically lasting between 30 and 60 seconds; i.e., similar to the period 
of the major wave groups. 

Changes in the wave field and the wind-wave relationships throu^i 
an ’event 1 are illustrated for the burst of segment 8, run 157, in Figs. 

25 and 26. Fig. 25 shows the coherence of U and rj from before, during, 
and after the event (segments 6-10); Fig. 26 shows the momentum spectra 
from the same period. Before the event is high and = 'IT for the 

low frequency components; is insignificant for most higher frequencies. 

During the event increases dramatically in the * 1 . 25-3 . 25 hz range. © u *j 
shows a continuous ’slip 1 through. this range; this is probably accounted 
for by the averaging time being much larger than the high frequency insta- 
bilities so that both generation and detritus is represented, and by the 
deceleration of the gust parcel during generation. Immediately after the 
event Torj has decreased for frequencies >1 hz;*it is still lower in segment 
10, except for the local peak at 1.25 hz where the waves lead the wind (an 
opposite phase relationship from that of segment 9). Before the event the 
wave momentum is strongly concentrated in the local waves (0.32 hz) and 
swell (0.15 hz). During the burst the high frequency components, out to 
over 3 hz, carry a much larger portion of the momentum than at any other 
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Figure 23 
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Figure 24. 
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time. Immediately after the event the very high frequency components have 
disappeared, and the 0.45 hz component has been strongly augmented. This 
peak shifts to a still lower frequency CO. 40 hz) in the following segment, 
and eventually nearly all the momentum is again concentrated in the local 
waves and swell. 

Fig. 20 and Fig. 21 show that the event, on the 30 second averaging 
basis, was notable by anomalously high values of U^q, wind shear,; wave height, 
current maximum, E, E(30 sec), and M. T increased slightly during the 
event, and significantly 60 seconds later. 'C ^ is not excessively large, 
but the profile for segment 8 had a low correlation with the model 
( P - 0.91) so that the computed value of ~C is unrealistic. The large 
departure of profiles from the logarithmic mean, and the development of 
unstable shear layers, are characteristic of such intermittent events (as 
noted in the cat’s paw study and further shown in Figs. 28 and 29). 

Fig. 27 depicts the detailed (analog strip chart output) structure 
of such a burst type generation event from run 116. Isotachs (from profile 
measurements), wind velocity fluctuations at two heights (high pass filter- 
ed at 0.05 hz) , wave height, and a measure of the energy content in low and 
high frequency bands of wind and waves are depicted. Scales are relative. 

The major burst occurs after a rapid wind increase and is associated with 
the coupled occurance of low frequency wave groups in air and water. The 
generation of atmospheric turbulence (very high frequency wind fluctuations) 
occurs predominantly over the wave peaks in air and water, and these peaks 
occur simultaneously. Both low frequency and high frequency wind and wave 
energy are maximum during the burst, and some intermittent activity con- 
tinues after the major gust had passed. The low frequency wave components 
are slightly damped after the gust leading edge, when the burst is strong- 
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Figure 27. 
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Figure 29. 
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est. This figure clearly demonstrates a very strong coupling between 
media and between widely separated frequency components. 

Figs. 28 and 29 are similar diagrams from the first and third E peaks 
(see Fig. 20) of run 159 (the scale is comparable in Figs. 27 and 29; 
all others are relative). These figures again show that intermittent, high 
growth rate wave and turbulence generation occurs during the coupled 
occurrance of large, low frequency groups in wind and surface wave fields. 
The event is preceeded by profile distortions, a wind decrease then rapid 
increase. In Fig. 28 the low frequency components are again seen to be 
slightly damped during the intense high frequency instabilities. Large, 
low frequency components of the wind and wave fields that occur indepen- 
dently do not lead to this intense coupling and momentum transfer to high 
frequency components; it is the strong nonlinear coupling between the 
media which results in the generation of high frequency instabilities and 
enhanced coupling of scales. 

The strong shear layers in the atmospheric profiles are of very 
short duration. They are almost invariable associated with large wave 
groups and high frequency wave energy spikes, and then quickly revert to 
logarithmic form as the mean wind energy is drained and the atmospheric 
dynamics recover from the local instabilities. The current profiles show 
the effects of the momentum of the large wave groups; they have little 
direct relationship to the instantaneous wind profiles. 

The coherence and phase diagrams for the last event of run 159 are 
shown in Figs. 30 and 31. Most of the large energy peak for this time 
period was due to the presence of an exceedingly large wave group, and 
the high frequency generation was less intense than in the similar set from 
run 157. Nevertheless, the sequence of events associated with the instabil- 
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